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A PREDICTIVE PILOT MODEL FOR 
STOL AIRCRAFT LANDING 

by 

David L. Kleinman and William R. Killingsworth 
Systems Control, Inc. 


SUMMARY 


The optimal control approach to pilot vehicle systems analysis provides the 
framework for modeling pilot performance during STOL flare and landing. The 
model includes both the terminal time aspects and the short-term open/closed 
interrelationship for the landing task, and pilot reaction to ground effect. 
Model output predictions include the probability densities of the touchdown time 
and velocity, flare path dispersions and pilot control inputs. The model is 
used to predict pilot landing performance for three STOL configurations, each 
having a different level of automatic control augmentation. Model predictions 
are compared with NASA flight simulator data. It is concluded that the model 
can be an effective design tool for studying analytically the effects of display 
modifications, different stability augmentation systems, and proposed changes 
in the landing area geometry. 
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UNIT CONVERSION FACTORS 


The numerical results presented in this report are in the U.S. system of 
units. The conversions between the U.S. units and the S.I. system of units is 
given below for those quantities appearing in this report. 
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I . INTRODUCTION 


Increasing traffic density and resistance to airport expansion have led to 
the consideration of short take-off and landing (STOL) aircraft. These aircraft 
are designed with special high lift devices that enable them to land at lower 
airspeeds with steeper glideslope angles than conventional aircraft. Less 
runway is needed and less ground noise is heard. Because of the steeper 
approach paths and smaller landing zones, the flare portion of the landing takes 
a much shorter time and the pilot is allowed little margine for error. Previous 
studies have indicated a high level of difficulty for STOL landing and wide 
dispersions in the touchdown data. Much of the difficulty is caused by adverse 
ground effects to which the pilot must adapt very quickly. 

This study addresses pilot performance for instrument approach and landing 
of a STOL aircraft, i.e., no external visual cues are available to the pilot, 
only the cockpit displays. Pilot models based on modern control theory are 
developed and used in conjunction with mathematical models of the STOL to pre- 
dict pilot performance in the longitudinal mode of the landing task. 

The STOL aircraft design investigated in this effort was a four engine 
subsonic jet transport with high-bypass-ratio turbofan engines and a high wing 
equipped with an external-flow jet flap. The engines were pod-mounted under 
the wing in such a manner that the jet exhaust impinged directly on the trailing- 
edge flap system. This arrangement provided the high lift required for short 
field operations. 

This basic aircraft possessed several inherent stability and control prob- 
lems. In the approach configuration, the equilibrium trim pitch attitude was 
-4.1°. Pilots disliked this condition due to the large change required in 
pitch attitude just prior to touchdown. Average pilot rating assigned to the 
longitudinal handling qualities of the basic aircraft was 6.5^, with the 
major objections being poor airspeed control, sluggish initial pitch response, 
large pitch attitude excursions associated with changes in thrust and flaps, 
a phugoid of short period that cuased pilot induced oscillations, and low 
apparent pitch damping. 



In attempting to overcome these control difficulties, several longitudinal 
augmentation schemes were investigated by NASA-LRC. One augmentation scheme was 
an autospeed system that drove the third segment flap to maintain a desired 

l - 

airspeed. The autospeed system accomplished three objectives. It eliminated 

• ■ i'i. 

the’phugoid mode that gave rise to many of the basic longitudinal handling 
qualities, it provided good speed control, and it relieved the pilot of the 
speed control task, thus reducing pilot workload. More importantly, since the 

autospeed system moved the flaps to maintain speed, the airplane could be 

■ i • ' r » ■ J ‘. r . 1 ■'•.>' \ 

trimmed to a nose-up attitude on the approach path. However, the pilot still 
had to actively control pitch to maintain this attitude. 

A second augmentation scheme incorporated a pitch attitude command and 
hold system with autospeed. The attitude command system allowed the pilot to 
trim the aircraft to the pitch attitude required for touchdown early in the 
approach. Once trimmed, no subsequent pitch control was required since the 
command system automatically maintained the desired attitude. The longitudinal 
piloting task was thus reduced to altitude control with throttles. 

In landing the aircraft, the pilot has a standard set of instruments. An 
eight-ball was used to display pitch attitude, an altimeter and sink-rate 
instrument displayed vertical information, and an airspeed indicator provided 
velocity information. A light was automatically turned on at an altitude of 
53* to notify the pilot at flare initiation altitude.* 

The STOL iarcraft approached the runway on a 6° glideslope at a velocity 
of 75 knots. The primary longitudinal requirements for a satisfactory landing 
were: touchdown within a prescribed landing zone 450 ft long starting 250 ft 
from the runway threshold; and touchdown with an acceptable rate of sink, no 
greater than 3 ft/sec. Since flare initiation is a h^ = 53 ft touchdown 
occurs in approximately 5 sec. As a result, the pilot has little time to 
correct any errors. The difficulty of the landing task is compounded by adverse 
ground effects that arise from the deflected thrust. Specifically, a nose-down 
pitching moment, a loss in lift, and a decrease in drag occur as the ground is 


An additional "get ready" light came on -5 sec prior to flare initiation. 
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approached. 

In the analytic effort described herein, mathematical models were developed 
for the STOL aircraft and for the pilot. These models were then used to predict 
pilot performance in the STOL landing task for the basic aircraft, aircraft 
with autospeed and the fully augmented STOL. The pilot model was developed using 
techniques of modern control and estimation theory, and contains explicit repre- 
sentations of the pilot's inherent limitations, information processing and 
control behavior. The model extends existing results in man-machine systems 
analysis by considering human adaptation to ground effects as well as the 
terminal aspects and the open/closed loop interrelationships. 

The outline of the report is as follows. In Chapter II, linearized models 
for the three different STOL configurations are developed from NASA supplied 
wind tunnel data. In Chapter III, the modern control approach to human operator 
modeling is outlined, and the extensions necessary to treat the landing task 
are presented. Expressions for the touchdown statistics are derived from model 
covariance predictions. For the three cases studied, numerical values for the 
pilot model parameters are chosen a priori. In Chapter IV, the resulting model 
covariance predictions are compared with ensemble averages of actual pilot 
response data obtained on the Real Time Dynamic Simulation Facility at NASA 
Langley Research Center. The comparisons were favorable. Discrepancies between 
the model prediction and the data are discussed, and are used to indicate 
potential areas for model improvement. Our conclusions and suggestions for 
further work are contained in Chapter V. 
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II. STOL MODEL DEVELOPMENT 


In this chapter, mathematical representations are developed for the STOL 
aircraft during approach and landing. In particular, linearized dynamics are 
derived from available wind tunnel data. The linear models are needed in the 
subsequent development of the optimal control model for pilot response that 
includes the terminal control aspects of the landing task. 

2.1 BASIC EQUATIONS OF MOTION 

The equations of motion for the aircraft can be derived from Newton's 
laws of motion. 


£F = 


m = 



( 1 ) 

( 2 ) 


where the subscript I indicates the time change with respect to inertial space. 

Assuming: (a) the mass of the aircraft remains constant during any partic- 
ular dynamic analysis; (b) the body fixed axes OX and OZ lie in the plane of 
symmetry of the aircraft; (c) the aircraft is a rigid body; and (d) the earth 
is an inertial reference, Eqs. (1) and (2) may be expanded to yield 

EF x = m(U + WQ - VR) 

EF = m(V + UR - WP) (3) 

y 

EF z = m(W + VP - UQ) 
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and 


EM = PI - RI 
x x xz 


+ QR(I -I ) PQI 

z y xz 


EM 

y 


QI + PR(I -I ) + (P 2 -R 2 )I 
y x z' xz 



( 4 ) 


EM = RI - PI + PQ(I -I ) + QRI 
z z xz y x xz 


These six simultaneous nonlinear equations of motion completely describe 

the behavior of a rigid aircraft. In this form, a solution cari be obtained 

only by use of analog or digital computers. In most cases,' however, by use of 

proper assumptions, the equations can be broken dowh. into two sets of three , 

equations each and these linearized to obtain equations amenable to analytic 

solutions of sufficient accuracy. The six equations are first broken up into 

two sets of three simultaneous equations. To accomplish this, the aircraft is 

considered to be in straight arid level unaccelerated flight and then to be 

disturbed by deflection of the elevator. This deflection applies a pitching 

moment about the OY axis, causing a rotation about this axis which eventually 

causes a change in F and F , but does not cause a rolling or yawing moment 

X z 

or any change In F ; thus, P=R=V=0 and the EF , EM^, £M z equations may be 
eliminated. This leaves three equations describing the so-called longitudinal 
motion: .■ r . '> 


EF x = m(U + WQ) ■ 

' J 

EF z = m(W - UQ) 


(5) 


( 6 ) 


EM = QI 

y y 

i 


(7) 


The applied forces and moments acting upon the aircraft will next be considered, 
and then the equilibrium longitudinal flight conditions will be examined. 
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2.2 APPLIED FORCES AND MOMENTS 


To complete the mathematical representation, it is necessary to expand 
the applied forces and moments which are of an aerodynamic or gravitational 
origin. For example, the components of gravity along the X and Z axes are 
a single function of the angle 9 between the X axis and the horizontal, 

F , gravity = mg sin 9 

( 8 ) 

F , gravity = mg cos 9 


For the longitudinal analysis, the forces in the X direction and the Z 

direction may be considered to be functions of: (a) the coefficient of thrust, 

C ; (b) the angle of attack, a; and (c) the controls fit = tail deflection, 

r 1 81 

and fif^ = flap deflection from 60°. These forces may be expressed as: ’ 


EF = 


» v a s r 


x 2 


pv 2 s 

EF = 
z 2 


C x (C T ,a) + C x (C T ,a)St + C x (C T ,a)6f 3 

O fit fif^ 


C (C ,a) + C (C ,a)fit + C (C_,a)6f 

z o 1 z 5t 1 z 5f~ 5 

J 


- mg sin 9 


+ mg cos 9 


(9) 


( 10 ) 


where the coefficients ^xg t > etc., are determined empirically by wind 

tunnel testing. (In this testing, thrust variation was studied by blowing 
compressed air through the engines of the model. The coefficients are tabulated 
and/or presented graphically in Ref. [1].) The Y-moment equation describing 
pitching motion is given as 



pV 2 Sc 

a 

2 


) 


C (C a) + C (C ,a) fit + C (C T ,a)Sf 
mi m. i j 

o fit fif^ 


+ 


2 2 
P V a SC 


C m (C T’ a)Q + C m- (C T’ a) “ 

L q a 


(ID 
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where the C coefficients are obtained from Ref. [11. The differential 
m 

equations describing the longitudinal dynamics are ;then given as 


U = -g sin 9 - WQ + 


2 

pV s - 

a 


2m 


C + C St + C Sf„ 
, x Xo £ 3 

L O St Of^ 


( 12 ) 


W = g cos 0 + UQ + 


pv 2 s - 

a 


2m 


C + C St + C 6f_ 

Z Z r. Zp £ <3 

o ot 


(13) 


Q = 


pV 2 Sc 

a 


21 


y l o 


C + C St + C 6f. 

m m f m » _ 3 


ot 


<5f. 


2 2 

pV^Sc r 
a 


i' i. ■ 


41 


C Q + C a 
m m* 

L q a . 


(14) 


The equilibrium flight conditions that satisfy these equations will next be 
determined. 

2.3 EQUILIBRIUM FLIGHT CONDITIONS 

For steady flight conditions, the rates of change in the dynamic equations 
are zero, i.e., U=W=Q=0. For the basic aircraft, we shall consider that the 

horizontal tail is the only control used to maintain equilibrium flight condi- 
tions and that the flap is set at full 60° deflection, i.e., Sf^ = 0. There- 
fore, the equilibrium conditions are described by 


0 = 


pvjs r 

-g sin 0 + 


C + C St 
, x x~. 

L o St . 


(15) 


0 = g' cos 0 + 


pV 2 S 

a 

2m 


C + C St 

L Z o z St J 


(16) 


0 = 


pV Sc r 
a 1 

21 


C + C St 
m 

- o ot . 


(17) 
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2.3.1 Steady Level Flight 


For straight and level steady flight, the flight path angle y equals zero. 
Thus, since y is defined as 


y = 0 - a. 


(18) 


for level flight we have 0 = a. The X force equation is then 


pV a S 

+ ~5T c x 'V 01 ' + c x , 5t - 0 

O Ot 


(19) 


For straight and level flight prior to initiating the glide slope, V =75 knots 
and Thrust = 23,000 lbs. This yields a coefficient of thrust 


C = — y- = 1.71 

pV S 
a 


( 20 ) 


In solving for a, (C T = 1.71, a) is approximated as 

o 


'9C 


x 


C x (C = 1.71,0) = C x (C t = 1.71, a=0) (C T =i.71,a=0) (21) 

o o 


where from curves of C , we obtain 

x 

o 


C x (C = 1.7. ,a=0) - (-.25) C x (C T = 1.71,ot=0) - 0.0 


For a within a region near zero, C may also be neglected. 

x 6t 


( 22 ) 


* 

The notation 




3a 


is used. 
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Hence 


pV 2 S 

“o - 2iST C * (C T ■ 1 - 71 >“-°> 

o 


(23) 


or 


a = (.284) (-.25) = -.071 rad = -4.1 deg = 0 
o o 


(24) 


1 i 


To maintain constant pitch, it is necessary to introduce a constant longitudinal 
trim 6t. This is calculated using the equilibrium Y-moment equation. 


0 = 


pV 2 Sc 

a 

21 


C (C ,a) + C (C ,a)<$t 
m l m r l 

. o ot 


(25) 


This requires 


or 


C (c = 1.71, a = -4.1°) = C (C m = 1.71, a = -4.1°)6t 
m T o m~ T . 

o ot 


6t * - 1 ' 75 


(26) 


Hence, for straight and level flight at V =75 kts, T = 23,000, the aircraft 
is pitched down at -4.1° and the horizontal tail has a trim value of -1.75°. 


2.3.2 Steady Flight Down Glideslope 


On the glideslope, the velocity is held at 75 knots and the glideslope and 
flight path angle are -6°. The relationship between flight path angle and angle- 
of-attack is given by Eq. (18). Thus, 0 = a - .105. The X-force equation is 
then given as 


0 = -g(a - 0.105) 



C x (C T ,a) 
o 


+ C (C ,a) 6t 
x St 1 


(27) 
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( 28 ) 


again neglecting ((^,,0=0) and C (C ,cs=0), 

x a x 5t 


pV a 

“■ - 105 + 2Si C Xo < C t-“"» 


Now noting for a range of C^, about 1.0 that (C^,,cc=0) is approximately 
constant at -.25 yields ,° 


a = 0.034 rad = 1.94' 
o 


and 


0 = -4.1' 
o 


The vertical equation in equilibrium is given as 


(29) 


0 = g cos 0 + 


PV 


2 


2m 


C (C a ) + C (C a )6t 

Z I o z. t I o 

o ot 


(30) 


Neglecting the contribution of the tail yields 


c z ‘W ' - 4 - 1 ' 

o 


(31) 


This condition for C is satisfied at C_ =1.32 or T = 17,500. From the 

z T 

o o 

pitch equation, the longitudinal trim, is given by 


-C (C ,<x ) = C (C ,a )St 

IQ i O 01 r> J. O 

O O Ot O 


yielding 


r — . 08 , c 

6t = T^8 = - 1 


( 32 ) 
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These equilibrium points analytically determined agree very closely to those 

r 9 1 

obtained via NASA simulation of the nonlinear equations and relationships. 

2.3.3 Steady Pitched-Up Flight Down Glideslope 

i 

In NASA simulations,^ subject pilots objected to the nose-down attitude 
because of the relatively large change in pitch attitude required prior to 
landing. In order to fly the plane in a pitched-up attitude, less than 60° of 
flap deflection must be used. An equilibrium pitched-up condition with a <5f^ 
input is now calculated. ' 

In the approach condition, the thrust is set at 13,200 lbs and the desired 
attitude is 2° pitched up. Since the aircraft is approaching on a 6° glideslope, 
the angle-of-attack is therefore 8°. For equilibrium flight, the X force 
equation satisfies the relationship 


0 


sin 0 + 


2 

pV S 
a 

2mg 


C + C fit + C fif 
X o x 6t x 6f 3 


Using coefficients from NASA data, this equation becomes 


(33) 


. 0048fit + .016 fif = -.224 . 
The moment equation must satisfy 


0 = 


pV 2 Sc 

a 

2iy 


C + C fit + C 5f, 
in m (> m pi p 3 

O Ot Of^ 


(34) 


or 

-.093 fit + .0025 fif = .3 . 

Solution of- these simultaneous equations yields fif 3 = -13.1° and fit = -3.5°. 
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2.4 LINEAR MODEL DEVELOPMENT 


Since the general nonlinear longitudinal equations are not amenable, to 
analytic study, linear models of the aircraft are now developed. In the analysis 
to follow, the aircraft is always considered to be in equilibrium flight before 
a disturbance is introduced. The linearized model of the aircraft then describes 
the dynamic behavior of the aircraft about the equilibrium operating point. 

2.4.1 Linear Equations of Motion 

Previously, the components of the total instantaneous values of the linear 
and angular velocity resolved into the aircraft axes were designated as U, V, 

W, P, Q and R. Since these values include an equilibrium value and the change 
from the steady state, they may be expressed as 

U = U + u 
o 

W = W + w 
o 

Q - Q 0 + q 

where U q , W q , etc., are the equilibrium values and u, w, etc., are the 
changes in these values resulting from some disturbance. 

Similarly, the external forces and moments are written as the sum of equi- 
librium components and deviation from equilibrium, e.g. , 

EF = IF + ZAF . 

XX X 

o 



For the development of the linear model, a body axes system known as the 
"stability axes system" is introduced. This system is widely utilized in air- 
craft analysis and introduces several simplifications. For example, the OX 
axis could be aligned with the longitudinal axis of the aircraft; however, if 
it is originally aligned with the equilibrium direction of the velocity vector 
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of the aircraft, then = 0. This is illustrated in Figure 1. In Figure 1 

the axes X , Y , Z are earth reference axes; the X , Y , Z , are equilibrium 
e e e o o o 

aircraft axes. For any particular problem the aircraft axes, after being 
aligned with the X axis into the relative wind, remain fixed to the aircraft 
during the study of perturbations from that initial flight condition. Such 
a set of aircraft axes are referred to as "stability axes". The stability axes 
will be used in all the dynamic analysis to follow. As can be seen from 
Figure 1, 0 q and y^ are measured from the horizontal to the stability X 
axis. The angle y is often referred to as the "flight path angle" and is 
defined as the angle measured, in the vertical plane, between the horizontal 
and the velocity vector of the aircraft. By using stability axes, 0 q and 
are equal. The definition of attack is standard, that is, the angle between the 
velocity vector (or the relative wind) and the wind chord. As the change in 0, 

which is equal to 0, is caused by a rotation about the Y axis then q = 0. 

• • 

Under these conditions U = U + u, W = w. and, as U . is constant, U = u 

. o o 

and W = w. As the aircraft is initially in unaccelerated flight, Q q must be 
zero, thus Q = q. Making these substitutions, the force equations (5) and (6) 
become 

EAF^ = m(u + wq) (36) 

EAF^ = m(w - U Q q - uq) (37) 

By restricting the disturbances to small perturbations about the equilibrium 
condition, the product of the variations will be small in comparison with the 
variations and can be neglected, and the small angle assumptions can be made 
relative to the angles between the equilibrium and disturbed axes. This last 
assumption somewhat limits the applicability of the equations, but reduces them 
to linear equations. Thus, Eqs. (36) and (37) may be written as follows with 
the addition of the pitching moment equation 


EAF = m u 

(38) 

X 

EAF = m(w - U q) 
z o n 

(39) 

I AM = I q 

(40) 

y y 
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It is now necessary to include the applied forces and moments in the 
linearized framework. From Eq. (9) , the force in the X direction is given as 


pV 2 S " 

EF + EAF ) = -g sin 0 H 

x x / 2m 

o 


- { EF + EAF ) = -i 

m \ X o x / 


C x (C T ,a) + C x (C T ,a)6t 
o fit 


+ C (C a)6f 
«3 


(41) 


Thus, the force in the X direction may be considered to be a function of 
a (angle-of-attack) , T (thrust) , and 0. Linearization about the equilibrium 
operating point yields 


- (zf 

m V x 
x o 


pV 2 S ( _ 

+ EAF x J = -g sin 0 Q - (g cos 0 q ) + -jf- j C x (C T ,a Q ) 

loo 


+(2/V )C u+C a + C fiT+ 
a X o x a x 6t 


C (C ,a ) + C a 
x p T o x ~ 
fit o fit 

a 


+ C 6Tj (fit + fit) + 

Xt J ° 


C K„ (C T - a o> * “ 

fif, o 6f 

a 


+ C 5T I (fif H 

J e ) 


(42) 


where 


a 


*&) 


(°T ,a o ) 
o 


(43) 


x 


fit 


*&) 


a 


( C T > a Q ) » etc * 
o 


(44) 
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Then since the equilibrium condition is 


1 'a 

- ZF = -g sin 0 + ■— C «L , a ) = 0 (45) 

mx o 2m x T o 

o o o 


and neglecting second order terms, the force in the X direction is 


i EAF 
m x 


+ (pV S/m) C u (46) 

a X 

O 


= -(g cos 


e o> 


pv 2 s 

a 

2m 


Ca + C 6T + C St + C 


6f, 


a 


ot 


<5t 


Sf, 


where 



C x (C'j, ,ot^) , etc. 
St o 


The differential equation describing velocity variation in the X direction is 
therefore given as 


u + (g cos © o )0 - 


pv 2 s 

a 

2m 


C a 
x 

a 


pv 2 s 

a 

2m 


C ST + C St + C Sf „ 

X ST x 6t x 6f 3 3 


+ (pV S/m) C u 
a x 

o 


(47) 


Similarly, the linear differential equations describing the variations in the 
Z direction and the pitching motion may be derived as 


2 

pV S 

U a - U q + (g sin 0)9 

o o o 2m 


[v] • 


pv 3 s 

2m 


C ST + C St + C 5f 


'ST St 


'5f, 


+ (pV S/m) C u 
a z 

o 


(48) 
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and 


pV Sc 2 pV Sc 2 pV 2 Sc pV 2 Sc " 

q - — 77 — C q - C a - C a = ~— 

4Iy m 41 m* 21 m 21 

q y a y a y 


C 6 t 
m <ST 


+ C St 4- 6 f 

m St 


m 6f . 3 


( 49 ) 


It should be noted that the coefficients C , C , C , C , C , and C , 

x„ x t z z rm m in, 

a St a 6 t a St 

are obtained from slope measurements of curves of C , etc. Also note that 


C = C (C ,a ) 
m m T o 
q q o 


C = C (C ,a ) 
m a m a T o 0 


2.4.2 Short-Period Mode Approximation and Model Validation. 

1. Short-Period Mode Analysis : The characteristic modes for nearly all 

aircraft in most flight conditions have two oscillations: one of short period 
with relatively heavy damping, the other of long period with very light damping. 
The "short-period mode" consists primarily of variations in a and 0 with very 
little change in the- forward velocity. Therefore, to approximate the short- 
period mode, let u = 0 in the equations of motion and also, since forces in 
the X direction contribute mostly to changes in forward speed, the X equation 
is neglected. The equation may then be put into the following form (using Laplace 
variables for this segment of the analysis) 


/— s - C \ a(s) + - — s + sin 0 1 0(s) = C 6 T + C St 

\Sq Z a / Sq Sq °J Z 6T Z 6t 


+ C 




(50) 
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i. rs 

- -^7 C s - C a(s) + — s - C s 6(8) = C <$T + C 6t 
2U m* m — 2U m m nu 

a a Sqc q oT ot 


+ C 6f 

® ^ r; ^ 

6 f 3 


( 51 ) 


- PV a 

where q = — . Letting 0 q = 0, the determinant of the homogeneous equation 
may be put into the form 


s(D„s +D 1 s+D)=0 
2 i o 


(52) 


where 


D 2 = 


'4?V*) 

k Sqc/\ Sq/ 




(53) 


D = ~=r C C - — C 
° 2U m q z a Sq m a 


Dividing Eq. (52) by and writing in the standard quadratic form with Z, 

and a) , then 
n 


2?u) 

n 



and 


2 


u 

n 


D 


o 



(54) 


Explicitly we obtain 


0 ) 

n 


UpSc 

2 



(55) 
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and 




i / 2I y ' 

4 C m + C „- + “I C z 
\ q a me a/ 


me 


' C C 2mC 
m z m 

q a a 


pSc 


(56) 


2. Linearized Model Validation : The STOL transport aircraft has been 

simulated at NASA-Langley Research Center using the general nonlinear equations 
with the aerodynamic coefficients obtained via table look-up. The short-period 
dynamic response characteristics of the aircraft have been determined empirically 
by measuring the response of the simulated aircraft to simple inputs. Results 
for the basic aircraft at 75 knots and T = 23,000 are. 


to = 1.41 r = 0.71 . 
n 

These numbers provide an opportunity to evaluate the linear model developed 
herein by calculating io n and £ using the derived short-period expressions. 

From the aerodynamic data, we obtain 



C 

m 

q 



,a o ) 


-27.00 



-7.45 


C 

tn 

a 



(C.J, >01^). 2.1 

o 




m* 


= C 


a 


m* 


(C 


a 


,a ) 
T * o 
o 


- 11.0 


Substituting these values into Eqs. (55) and (56) yields 
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to 

n 


1.44 


? = 0.66 


These results agree very well with those values obtained experimentally with 
the nonlinear simulation and lend credence to the linearized model. 

2.4.3 Development of State Space STOL Models 


The linearized differential equation describing velocity variation in the 
X direction given in Eq. (47) may be rewritten as 


u' = (2qSC /mV )u' + (qS/mVa)a' - (g cos 9 /V ) 
x a o a 


+ (qS/mV ) 

cL 


C fiT + C fit + C fif„ 
x fiT x fit x 6f 3 3 


(57) 


the equation describing variation in angle-of-attack as 


a' = (2qS/mV ) C u' + (qS/mV )C a’ + q - (g sin 6 /V )0 
a z a z o a 

o a 


+ (qS/mV ) 

a 


C fiT + C fit + C fif_ 
. Z fiT Z fit Z fif 3 3 


(58) 


and the pitching moment equation as 


q + (qSc/I )C a’ 

y a 


(pV Sc /4I )C 
a y m 


a' 


a' + (pV Sc /4I )C q 
a y m 

q 


= (qSc/I ) 

y 


C fiT 4- C fit + C 
m fiT “fit m fif 3 



(59) 


where we have introduced the variables 

u' = u/V a' = w/V q = (pV 2 /2). 

& a. a 
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Substituting the expression for a' given by Eq. (48) into Eq. (59) yields 


q = (pV Sc /4I ) (2qS/mV )C u’ + (qSc/l )C . 

a y a z y m . 

J o ( 3 a 

+ (pV Sc 2 /4I )(qS/mV )C C ia + !(pV Sc 2 /4I )C 
K a y a m* z ( a 'ym 

J a a) ( J q 


+ (pV Sc 2 /4I )C (q - (pV Sc 2 /4I )(g sin .0 /V )C jo 
a y m a ) ( a y ° a m a) 


+ (qSc/I y ) 


C <$T + G 6t + C 6f 0 
iq r._- m ' in a ^ .3 

St St - Sf 3 


+ (pV Sc /4I ) (as/mV )C 
a y a m* 

a 


C ST + C St + C Sf 
Z ST z 5t Z Sf 3 3 J 


Now defining a state vector x = col[u’, a’, q, 0] and a control vector 
u = col [ ST, St, 6f 3 ] - Equations (57), (58) and (60) may be put into the 
matrix form 

x(t) = A x(t) + B u(t) 

where the matrix elements a(i,j), b(i,j) are 
a (1,1) = (2qSC /mV ) 

X 3 

O 

.... a(l,2) = (qS/mV ), ....... 

3 . 

a (1,3) = 0.0 

a(l,4) = -(g cos 0 /V ) 

O o 

a(2 , 1) = (2qS/mV ) .. 

Si 

a(2 ,2) = (qS/mV )C 

3 Z 

a 


(60) 


vector 


(61) 
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a(2,3) = 1.0 


a(2,4) = -(g sin 6 /V ) 

o a 

a(3,l) = (pV a Sc 2 /4I y )(2qS/mV a )C z 


a(3,2) - (qSc/I )C + (pV Sc 2 /4I )(qS/mV )C C 
y m a y a m z 

J rv J n 


a(3,3) = (PV Sc 2 /4I )C 4- (pV Sc 2 /4I )C 

y q y a 


a(3,4) = -(pV Sc 2 / 41 ) (g sin 8 /V )C . 

a y o a m* 

.< a 


a (4,1) =0.0 
a(4,2) = 0.0 
a (4, 3) = 1.0 
a(4,4) =0.0 

b(l,l) = (qS/mV )C , etc. 

a x St 

The remaining elements of the IS matrix are readily obtained from Equa- 
tions (57) - (60). To obtain models for various conditions, the aerodynamic 
coefficients are simply evaluated for the particular flight condition of 

interest. Stability derivatives such as C are calculated from slope 

m 

a 

measurement. 

To investigate approach and landing, it is necessary to include altitude 
variations within the model. Since the linear model describes variations about 
an equilibrium flight condition coming down the glideslope, the altitude 
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perturbation equation will describe altitude variations off the glideslope at 
a particular instant. Specifically, the sink rate is given as 

h = V sin y (62) 

and from the definition of the stability axis system, 

h = V sin (0 + 0 - a') (63) 

o 


Assuming small angles, the linear perturbation equation for vertical distance 
off the glideslope is given as 

5h = -Va' + V6 + u0 (64) 

o 

or dividing by V, 

Sh' = -a' + 0 + u'0 . (65) 

o 

One additional variable must be included to model the fact that changes in 
throttle do not lead to instantaneous changes in thrust. The dynamics of the 
engines are modeled as a first order lag in our linearized dynamic model of the 
STOL. Specifically, 

<5T = a„6T + a„6T (66) 

T T c 

where is the throttle commanded change in thrust level. Curves of 6T 

response to step changes in the <5T are given in Ref. [1]. From this data 

c _i 

an approximate thrust time constant a^ = .5 sec was selected. 

2.4.4 Wind Gusts and Ground Effect Disturbances 


Wind Gusts : The STOL aircraft is subject to vertical random turbulence 

during the approach phase. This is modeled by passing Gaussian white noise 
through a first order shaping filter, i.e.. 
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( 67 ) 


a + 

g 



where w is vertical wind gust velocity, a' is the increment in angle-of- 
g g 

attack created by the wind gust and where E(n(t)n(x)} = 6(t-t). L is taken 

w 

to be a nominal value of 500 and is 2, 4, and 6 for light, moderate, and 

heavy turbulence, respectively. (These numbers are chosen to correspond with 

NASA simulation parameters and, as in the simulation, the vertical gusts are 

reduced to zero between 100 and 50 feet altitude.) In studying the flare and 

landing portion of the flight, the turbulence may be neglected, since flare 

begins at h„ - 41 ft. (h - 53 ft.).* 

£g eg 

Ground Effect : As the externally blown flap STOL approaches the ground, 

adverse ground effects arise. Specifically, a nose down pitching moment, a loss 

in lift, and a decrease in drag are created. The effects of ground proximity 

on C and C begin at an altitude of approximately 60 ft and increase in 
o m 

severity, almost linearly with (60 - h) , to touchdown. The effects of ground 
proximity on C begins at a lower altitude and increases rapidly in an almost 

Xj 

linear manner. The incremental changes in pitching moment, lift, and drag 
coefficients due to ground effect are presented in Figure 2. 

The dynamic changes in the aircraft states due to the ground effect input 
may be expressed as 


The turbulence, by acting during approach, does effect the initial conditions 
at the flare altitude. This indirect influence of the turbulence is considered 
setting up conditions at flare initiation. 
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a 


ge 


C (h(t) ) 
mV D ge 




C (h(t)) 

. I M ge 

L y 


( 69 ) 


However, assuming a nominal value for C Q (h) =.-0.15 and a. nominal flare time 
(-5 sec.), the change in airspeed due to ground effect is a small increase of 
approximately two and a half feet per second. This small effect is neglected 
in the subsequent modeling. 


The lift and moment coefficients may be approximated as 


1 

< 0.0 


4.1 < h 


C L (h(t)) = 

[ (.0045) (h(t) 

- 41.) 

28. < h < 41. 

(70) 


l (.014) (h(t) - 

32.) 

0 £ h £ 28. 


! 

C M (h(t)) = ' 

0.0 


60. < h 

(71) 

(.004) (h(t) - 

• . s 

60.) 

0 £ h < 60. 



To include the ground effect within the dynamic model, two states are defined as 


x x (t) = C L (t) 
X 2 (t) = °M (t) 


0 £ h £ 41. 
0 < h < 60. 


( 72 ) 


Thus, the perturbations due to the ground effect may be expressed as 
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0 < h < 41. 


a' 


r qs x i 


mV 1 



o 

• 

_q_ 


ISc 

I 2 

ge 

L y J 


( 73 ) 


0 < h < 60. 


and the variations in and given as 


V 


m(t)h(t) 

0 £ h £ 41. 


= 


(74) 

A_ 


.004 h(t) 

0 £ h £ 60. 


where m(t) = .0045 for 28 < h < 41, .014 for 0 < h < 28. The total sink 

rate is 

h(t) = V y + V (-a' + 0 + u'0 ) (75) 

o o o o 

Equation (74) may be most conveniently modeled as a deterministic time varying 
input given as 


x x (t) 


z^a' (t) , 9(t), u*(t)) 


z- L (t) 

x 2 (t)_ 


z 2 (a’(t), 0(t), u*(t)) 


z 2 (t) 


Equations (76) and (74) may then be used to include ground effects as an external 
input disturbance in the STOL model. 


2.4.5 Dynamical Models of the STOL t 

Basic Aircraft : An overall dynamic model of the basic STOL aircraft 

during flare may be obtained by combining the equations describing the aircraft, 
and the ground effects. This model is given as 
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x=A}£+]5u + F^ z^(t) (77) 

where x = col[C C , u' , a', q, 0, <5h’, 6T] , u = col[6T , <St] 

lj M C 
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a il 

a l2 
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b ll 
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1 
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a 31 
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b 31 
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0 

0 

0 

0 

0 

0 

0 

0 

" a T 


a T 

0 


and 


F* 


10 0 0 
0 10 0 


0 0 0 0 
0 0 0 0 


(78) 


where 



z 1 (t) 

z x (t) 

z 2 (t) 


. 0045 [-V Y + V (-a' + 0 + u'0 )] 
o o o o 

. 014 [-V y + V (-a' + 0 + u'0 )] 
o o o o 

. 004 [-V y + V (-a' + 0 + u'0 )] 
o o o o 


28 <h < 41. 
0 < h < 28. 
0 < h < 60. 


(79) 
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The coefficients are evaluated for the aircraft in the basic pitched-down 
equilibrium condition. 

Autospeed : Through simulation experiments, the basic STOL aircraft repre- 

sented by Eq. (78) was found to possess several major longitudinal deficiencies: 

• sluggish initial pitch response; 

• low apparent pitch damping; 

• large pitch excursions associated with changes in thrust, flaps, and 
spoilers; and 

• a phugoid with an unusually short period. 

Due to these characteristics, there is poor pilot control of pitch attitude 
and, hence, poor control of airspeed. As a result of these difficulties, several 
augmentation schemes have been developed to improve the stability and control 
characteristics of the aircraft. 

The first augmentation system to be considered was an autospeed system 
that maintained the desired airspeed by driving the third segment flap. To 
neutralize the lift increment resulting from flap deflection, an interconnect 
to the symmetric spoilers was provided, thus effectively decoupling the forward 
and vertical modes. The autospeed system accomplished three objectives: (1) it 
eliminated the phugoid mode which was the source of much of the basic longitu- 
dinal handling difficulties, (2) it provided good speed control, and (3) it 
relieved the pilot of the speed control task and, hence, considerably reduced 
pilot workload. In addition, since the autospeed system moves the flaps to 
maintain speed, the pilot is able to trim the STOL in a nose-up attitude (as 
described in Section 2.3.3). It should be noted, however, that the pilot must 
still control pitch. 

A dynamic model describing the state variables of the STOL aircraft with 
autospeed is developed by removing forward velocity variation from the model 
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presented by (78). This is a justifiable modeling approach because of the 
decoupling nature of the autospeed system. This model is given in the form of 



(80) 


where the coefficients are evaluated for the pitched up equilibrium position. 

Autospeed plus Pitch Command and Hold : Due to continued difficulty in 

controlling pitch, a pitch attitude command system was incorporated with the 
autospeed. This system allowed the pilot to trim the aircraft to the pitch 
attitude required for touchdown early in the approach. No subsequent pitch 
changes are necessary since the command system automatically maintains the 
attitude throughout the remainder of the approach and landing. This configura- 
tion of the aircraft is modeled by removing pitch rate, pitch, and horizontal 
tail from Equation (80) . This yields the following model for the fully augmented 
STOL aircraft 
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where the coefficients are evaluated for the pitched-up equilibrium condition. 
In this configuration, the pilot's only task is to control altitude and sink 
rate with thrust. 

In summary, state space models in the form 

x(t) = A x(t) + B u(t) + F z^(t) (82) 

have been developed for three different STOL configurations: 

1. Basic Aircraft, 

2. Autospeed, and 

3. Autospeed plus pitch command and hold. 

Numerical values for the A and 15 matrices are given in Appendix B. In the 
next chapter, a model for the pilot will be developed for the flare task that 
can be used to predict pilot response. 
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III. PILOT MODEL DEVELOPMENT 


A model for the pilot is developed, based on the optimal control model for 
r 2—5 1 

human response, 1 that includes the terminal control aspects of the landing 
task and pilot adaptation to the ground effect. The treatment of the finite 
time aspects and open/closed loop interrelationship is an important extension 
of existing results in man-machine systems analysis. 

3.1 PILOT MODELING - GENERAL BACKGROUND 

The pilot modeling techniques used in this study are rooted in optimal 
control and estimation theory. They are based on the assumption that the well- 
trained , highly motivated pilot behaves in an optimal manner subject to his 
inherent limitations and constraints. This modeling approach is capable of 
treating multivariable, and time-varying systems within a single conceptual 
framework using state-space and time-domain techniques. It is, therefore, well- 
suited to the pilot-STOL system described in the previous sections. The 
generalized modeling framework facilitates extensions to cover the terminal 
control aspects of pilot response. 

The basic operator modeling techniques are documented extensively in the 
r 2—5 1 

literature; 1 therefore, only the salient features are described below in a 
general context. The model extensions necessary to treat the terminal time 
aspects of the landing task are developed in the following sections. 

3.1.1 System Description 

The structure of the optimal control model of human response is shown in 

Figure 3. It is assumed that the system dynamics (which also include any per- 

'fc 

tinent noise shaping filters) are described by the linear time-invariant 
equations 


The time-varying case is discussed in Ref. [5]. For the STOL problem studied 
here, the linearized system dynamics are time-invariant. 
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( 83 ) 


x(t) = A x(t) + B u(t) + £ w(t) + F z(t) 
x(0) = given 

where the n-vector x(0 represents the system state, u(t) = coltu^,^* . • .u^] 
are the human's control inputs, and where w(t) represents random input forcing 
functions, e.g., wind gusts. w(t) is assumed to be a zero-mean, Gaussian 
white noise process (possibly nonstationary) with covariance 

E{w(t) w '(a)} = W(t) S(t-a) (84) 

The term z^(t) in Eq. (83) is used to generate input forcing functions 
that are not stochastic in nature, but are more aptly regarded as "deterministic" 
inputs. In the modeling of pilot response, it is assumed that z(t) is the 
time rate of change of inputs such as wind shear, ground effect, etc. For example, 
the scalar equation 

x x (t) = z(t) (85) 

can be used to generate a "deterministic" input x^(t). We assume that the pilot 
can continuously estimate the quantity x^(t) = z(t) from displayed information, 
but does not so estimate z(t). Thus, x.. (t) must be modeled as a state variable 
as in Eq. (85). 

The human observes a set of outputs £(t) = col[y^,y 2 » . • . ,y m l that is 
related linearly to the system state and control, 

^(t) = C x(t) + D _u ( t) (86) 

where the matrices J3 and D can be time varying to model gains on display 
variables that change with time or altitude. The usual assumption in the model 

A generalization to Eq. (85) is x^(t) = X 2 (t), Xj(t) = z(t). This would be 
used of x^(t) appeared explicitly in the system equations (83) or contributed 
directly to a displayed signal. 
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is that if a quantity y^ is explicitly displayed to the man, he also derives 
the rate of change y^, Thus, ^(t) contains both position and velocity infor- 
mation of a displayed signal, but no higher derivative information. 

3.1.2 Human Limitations 


The human has inherent limitations of time-delay and perceptual noise 
(i.e., remnant). These quantities are associated with the observational process 
in the man model, so that the human is assumed to perceive ^(t), a delayed, 
noisy replica of y(t). Thus, 

y (t) = £(t-x) + VyCt-x) . 

The time-delay x is nominally x = . 2 + . 05 sec. The "observation" noises are 
white, independent, and have covariance 

E{v .(t) v .(a)} = V . (t) 6(t-a) i=l,2,...,m (87) 

yi yi yi 

When directly viewing y^(t), the 
scale with the variance of y^(t), 

V . (t) = p .E(y?(t)} (88) 

yi yi l 

For full attention on a single display indicator, the noise/signal ratios p ^ 
on position and rate typically have a value of .OItt, i.e., v ^(t) has a 
-20dB normalized power density level. When there are K > 1 display indicators, 
the human must allocate his attention among the various displays. Let n k 
denote the attentional allocation to display indicator k. Then, neglecting 
the time spent in inter instrument scanning, we have 


associated covariance V . is assumed to 

yi 


i.e. 


K 

2 n k = i, 0 < n k < i (89) 
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If displayed variable y^(t) is obtained from indicator k, the effect of 
attentional allocation is to modify the noise/signal ratio according to 





( 90 ) 


where - .Olir is the noise/signal ratio that corresponds to full attention 

on indicator k. Methods for determining within the optimal control model- 

ing context are discussed in Refs. [2] and [5]. However, these methods are 
difficult to apply, and are generally applicable for steady-state situations 
only. In the present STOL effort, we do not attempt to solve the display atten- 
tional allocation problem in the time-varying or finite time case. However, in 
order to include some effects of attentional allocation, albeit in a crude 
manner, we set = 1/K. Thus, we assume an equal division of pilot attention 
among the primary instruments needed for control. 

In addition to time-delay and perceptual noise, we include perceptual/ 
indifference thresholds on displayed information. Clearly, a pilot will not 
react to changes in a displayed variable if the indicator motion is smaller than 
his observational thresholds. Thus, in perceiving quantity y^(t), we write 

y (t) = f . (y. (t-x)) + v (t-x) (91) 

pi li yi 


where the threshold nonlinearity 


£ i<-> 


is 



x 5^ a 
-a _< x a 
x < -a 


(92) 


The nonlinearity f(«) is replaced in the pilot model by its equivalent 
"Random Input Describing Function". This statistical linearization gives 


f (x) - f • x 


(93) 
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where f is a gain that varies with the mean and standard deviation of x(t). 

X 

Expressions for f^ are given in [4], Values for the thresholds a^^ depend, 
in part, on the physiological limitations of the eye in sensing motion, as well 
as on the specific instrument markings used. 

3.1.3 Task Definition 


It is assumed that the human's control task is adequately reflected in the 
choice of a control u(*) that minimizes the cost functional 


J(u) = lim E — 
T-w» 1 


T 

[j_' (t) QyX(t) + u' (t)Q u u(t) + u' (t)Q^u(t) ]dt 


(94) 


conditioned on the perceived information y^(*)* The first term in J(u) is a 
generalized mean-squared error criterion where ^ depends on the task specifics. 
The control rate term is used to account for the human' s limitation on the rate 
of control motion, and introduces "neuro-motor" dynamics in the man model. 


The selection of the weightings = diag [q ], = diag [q ui ] and 

= diag [q^] in J(u) is a non-trivial step in applying the man-model. In 
any specific situation, these parameters are dependent on both the human's 
objective task requirements and his subjective mode of behavior. One useful 
method for selecting reasonable a priori estimates for q ^ and q^ is by 
associating these quantities with allowable deviations in the system variables. 
Thus, we let 



1 

2 



1 

y . 

i,max 


9 

^ui 

u. 

i,max 


(95) 


where y. is the maximum desired, or allowable, value of y u. is 

- / i,max i i,max 

the maximum control deflection. This method of choosing weightings has intuitive 
appeal. First, maximum or limiting values of system quantities are often easy 
to specify, or elicit by pilot questionnaire. Second, the contribution of each 
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is to its maximum 


term q .y7(t) to the total cost depends on how close y. 

yx x x 

value. Because of the normalization (95), each y^(t) will not contribute 


significantly to J(u) provided 

2 


lYil 


— ^i,max 


When 


|y i 


x,max' 


the 


contribution to J(ij) of q y (t) increases rapidly with y^(t). The analogy 
to manual control is that there will be little concern over minimizing y^(t) 
if this quantity is well within allowable limits. 


One method that can be used for selecting the weightings q*_^ on the 

human's control rate is similar to that for q . , i.e., let 

yx' ’ 


q* . = 
ux 


u. 


x,max 


(96) 


Here, u. is the maximum rate that a human can (or will) input a control 

x,max l 

u.. However, we choose to associate the weightings q* with the human limita- 

tions directly. Weighting u in the cost functional results in first order 

lags being introduced in the man model. There is a first order lag time constant 

T n . that corresponds to each control rate weighting q^V ; the smaller one sets 

q^V , the smaller is the resulting T n ^. In the model, the lags T ni are 

associated with the man's "neuromotor" dynamics, where past modeling efforts 

[ 2 ] 

show typically T - .1 sec. Thus, the weightings q*. are adjusted 

1 * ul 
iteratively until each - .1 sec. ^ 

3.1.4 The Pilot Model 


The "human's" control input (i.e., the control that minimizes J(u)) is 
generated by the feedback law 


T u(t) + u(t) = -L x(t) + v (t) = u (t) 
— n — — — u — c 


+ -u(t) 


(97) 


i 

If the resulting q^v weighting is such that l/\/q^ is much greater'than the 
physical rate at which one can move control u^, then Eq. (96) must be used. 
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where x(t) is the best estimate of the system state x (t) . v^(t) is 3 white 

"motor-noise" that represents a human's imprecise knowledge of generdtdd control 
inputs and has covariance 

v u (t)v^ (a)J - -V (t) 5(t-cr) i=l,2,..;,r (98) 

i i ’ U i i 

where a good approximation to the covariance V (t) is 

u i 

V (t) = p • E{u 2 (t)} i=l,2, . . . ,r (99) 

u. u. 1 

r i 

The motor noise/signal ratio p - .003ir, i.e., v (t) has approximately 

i i 

-25dB normalized power density level. 



The feedback gains L and the r x r matrix T are given by 


.-1 


where 


T = P 00 Q- 
— n —22 


- -22 -12 


( 100 ) 


P = 


p 1 

-11 i 
1 

— 1 2 * 

-J- 


p ? J 

p„ 

-12 I 

-22 


satisfies the equation 


£ n A + A'P U + C'SyC - p^- 1 P i2 - 0 


P n l + A'P 12 + C' 2y £ - P^' 1 l 22 ' 0 


(101) 


-1 


r 12 B + B« P 12 + D'S^D + ^ - P 22 Q^ P 22 - o 
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The, estimate x(t) is obtained . from the cascade combination of a Kalman 
filter and predictor that compensate optimally for the human's observation noise 
and time-delay, respectively. The filter is defined in terms of the augmented 
, state vector co1 [x(t), u(t)]. The filter generates £(t) = Ef)£.(t-x) |y_(o) , 

<J _< t} at time t x of the past state x( t-T ) and control u(t-x) from 


~ £(t) = A £(t) + B u^t-x) + l C , V y 1 [£ p (t) - C £(t)] 


( 102 ) 


where 





C I 


(103) 


and 


V 

-y 


is given by 


S/t) - <iiag ( v yl <t)« yl l 


(104) 


The matrix E!(t) , t t satisfies the time-varying equation 


I(t) = A E(t) + Z(t)A’ + W(t) - 




C E(t) 


(105) 


where 


W(t) = 


EWE’ + FZF' 



0 


T 1 V T' 1 
— n — u — n J 


(106) 


and Z(t) = diag [z^(t), z^(t),...,z^ (t)] (n = no. of deterministic inputs). 

i z n z 

z 

The initial conditions on Eqs. (102) - (105) are 
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(107) 


2.( T ) = io , I(T) = 4 

where is the man's a priori estimate, of the mean state E{ >^( t ) } at t=0. 

£ is the man's initial uncertainty in this estimate. 


The predictor generates the best estimate of x^) denoted by x^) from 
the estimator output jp(t) according to 

£(t) = A i(t) + B u,(t) + e^ T Z C'r 1 [^(t) - C £(t)] 

(108) 

£(t) = e^ T £(t) 


The above equations (83) - (108) define completely the input-output model 
of the pilot-vehicle system. Note that the model is time-varying since the 
matrix Z varies with time. These equations can be used to predict pilot 
response in closed-loop tracking tasks. The model inputs include the vehicle 
dynamics, cost functional weightings and human limitation parameters. The model 
outputs include statistical predictions of system performance. In steady-state 
situations, one may also compute input-output transfer functions and power 

r 2 1 

density spectra. 1 

3.2 COVARIANCE EXPRESSIONS AND STATISTICS 

Analyzing the behavior of a time-varying system is best accomplished using 
covariance propagation methods. This results in predictions of both the mean 
system response and the standard deviation in this response. The mean response 
is the system response to a specific input z^(t) in Eq. (83), and represents 
the result one would expect to find by (ensemble) averaging the results of many 
trials, each with the same x(t)* The standard deviation results from the 
randomness v^(t) and v^(t) that the human injects into the loop, as well as 
from any gust inputs w(t). 


If x(b) 0 and W(t) W = constant, then the man-model reaches a time- 
invariant steady-state. 
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In this section, we present equations for the propagation of ^(t) = mean 

state and the covariance X(t) = E{[^(t) - " X(t)]'). These 

expressions give us the capability of analyzing the statistics of the closed- 

2 2 

loop system response, and are needed to compute E{y^(t)} and E(u^(t)l in 
Eqs. (88) and (99). * 


The required results are most easily derived using the equations for the 
estimation error (see Ref. [4]) 


e^t) = X(t-T) - £ (t) 


(109) 


and the predicted estimate ^(t). These quantities satisfy 


€: (t) = A f e, (t) + K(t)v (t-r) + £ £(t-x) + v (t) 
_L I X y j “U 


( 110 ) 


i(t) = A±(t) +K(t) [C e x (t) +W y (t-T)) 


( 111 ) 


where 


Af = A - K = closed-loop filter matrix 


A = A - _B L = closed-loop control matrix 


K = E C'V ^ = estimator gain at time t 

y 


K = e AT K 


The total state 2 £(t) at time t is given by 


At 


X(t) = x(t) + e- e^t) + e 2 (t) 


(112) 


where e^Ct) t * ie P re< ^ act: *- on error. 
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Expressions for the covariances 

X(t) = E{[£(t) - £(t)] * [£(t) - 

E (t) = E{ [e . (t) - e . (t) ] • [e . (t) - e.(t)]'} ; i=l,2 
1 -l -l -x -i (119) 

X(t) = E{[x(t) - x(t)] • tx(t) -xCt)]'^ 

= X(t) + e^ T M(t) + M' (t) e^' x + E^t) e^' T + E 2 (t) 


44 



where M(t) is the cross covariance between _X.(t) and e^(t), are derived by 
subtracting Eqs. (115) - (117) from Eqs. (110) - (111) and (113) and taking the 
autocovariance- The results are a coupled set of linear matrix equations 


E (t) = ApE, (t) + E 1 (t)A’ + K(t) V (t-T) K’ (t) + B V B’ 

—1 — r— 1 —1 — t — —y — ir - 

M(t) = -yi(t) + M( t ) + K(t)C [E (t) - Z(t)] (120) 

X(t) = A X(t) + X(t)A’ + K(t)C M* (t) +M(t)C'K' (t) + K(t)V (t-x)K' (t) 

^ c y 


E 2 (t) = 




v.W)!' ^ (t - 0) do 


( 121 ) 


t-T 


where ^(t) satisfies Eq. (105). The initial conditions on X, E^, M, are 
taken as 

X(t) = M(t) = 0, E^x) = X(0) = given 
Finally, the mean and covariance of ^(*) at time t are simply 


£(t) = £ x(t) 

(122) 

Y(t) = cov[£(t)] = c X(t)C' 

(123) 


Efficient computer programs have been developed for integrating the mean and 
covariance equations (114) - (123) presented above. 

The above expressions for the signal means and covariances are of special 
interest in predicting the probability densities of system variables. Since 
the white noise random inputs w(t), v^(t) and v^(t) are assumed to be Gaussian, 
and since the closed-loop system is linearized, the system states x^,...,x n 
are Gaussian random variables. The x^ are non-stationary when the Kalman 
filter gains K in Eq. (110) are functions of time. 
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For a Gaussian variable x, the mean m and variance 
statistics for defining the probability density of x. 


are sufficient 


p(x) = — — exp 

/2tt a 


(x-m) 2 ) 
2 a 2 ) 


(124) 


Thus, for any state variable x.(t), with mean x.(t) given by Eq. (114), and 

2 1 1 
variance 0^(t) = X^(t) given by Eq. (119), 


p(x i (t)) 


1 

2tt o . (t) 
i 


exp 


[x. - X.(t)] 2 
2a 2 (t) 


(125) 


is the probability density function of x^ at time t. Similar expressions hold 
for p(y jL (t)) and p(u i (t)). 


Equation (125) is the univariate density of x^(t). Since the x^(t) are 

correlated, the multivariate density of x = (x 1 x ) is given by the 

r 61 i / n 

general expression, 1 1 

p (x- ,x, x ) = expj- -^[x-x(t)] 'X 1 [x-x(t)]j (126) 

t 1 2 n 2ir n |x(t)| ' 2 - -- ) 


where X is the state covariance matrix and |x| = det X. The multivariate 
distribution of = [y^ > 5 ^ 2 » • • • >y m 3 can also be written using the expressions 
for y^(t) and Y in Eqs. (122) - (123). 

3.3 TERMINAL CONTROL 

The above section described the pilot model that has been developed to 
treat regulation (i.e., tracking) tasks. Modeling pilot response during flare 
requires that we treat the terminal control, or finite time aspects, of the 
task — namely, to land the aircraft at some intended point on the runway, with 
a desired sink rate and attitude. 
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For the aircraft system defined by 

x(t) = A x(t) + B u(t) + E w(t) + F z(t) 


j(t) = C x(t) + D u(t) 


(127) 


with 

x(t=0) = x(0) = given 

it is assumed that a general set of terminal conditions 

H x (T ) + c = 0 ' ■ (128) 

d — 


must be satisfied where 



is the intended touchdown time. 


In order to extend the human operator model to include the terminal control 

task, it is assumed that the human generates a control input u (t) such that 
■ • -• , ~° o 

the response x (t) of the unforced system. 


x (t) = A x (t) + B u (t), x (0) = 0 
— o o o — o — 

^(t) = C 2^(1) + D Uq ( t ) 

meets the terminal conditions 

H x (T ) + c = 0 
o d — 


(129) 


(130) 


Clearly, many controls can accomplish this transfer. However, we assume that 
the pilot's control is the one generated with least control effort. Thus, we 
require u (t). to minimize 


r 

•/ 


u' (t) Q*u (t) dt 


(131) 
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where the cost functional weightings q- . on control rate are chosen as discussed 

in Section 3.1.3. The requirement that u^ minimize J q was motivated by the 

fact that a pilot tends to control an aircraft smoothly, minimizing any unneces- 

* 

sary rapid control motions. 

The control u (t) that meets the terminal constraint (130) while minimizing 
— o 

the "smoothness" criterion (131) is derived in Appendix A. u^t) an °P en ~ 
loop (i.e., feedforward) control given by 


u (t) = - Q* 1 B h(t) 


(132) 


where the time function h(t) is generated by 


h(t) = - A' h(t) 


(133) 


with the boundary condition 


h(T ) = H' [H W(0,T,)H' ] _1 c 

— d o — cr- d — o — 


(134) 


where 


W(0 


,T ,) = f e^° (T d 0) B Q* 1 B' e^o (T d 0) da 
’ d I — o^Ha — o 


(135) 


The matrices A , B , H are 
— o — o — o 


A more general cost functional 



0 


could also be considered. In the present effort, Eq. (131) seemed most reason- 
able in view of general pilot technique. 
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(136) 


The control u^(t) and the resulting state trajectory x^t) and displayed 
outputs ZqCi) meet the terminal condition (128) only when there are no external 
input forcing functions, wind gusts, or pilot randomness entering Eq. (127). 

The fact that there are such external disturbances that act on the system give 
rise to deviations between the actual x(t)> as generated by Eq. (127), and the 
nominal x (t) . Let 


6 x = x(t) - x (t) 
— o 

5 X = x( fc ) " ^(t) 


(137) 


denote these deviations. Also, let the pilot's control input be 
jj(t) = u (t) + 6 u(t) 


(138) 


where S u(t) is the pilot's corrective control action. The pilot generates 
u(t) to keep variations j5 x(t) from the nominal flight path "small". 


From equations (127) - (129), x(t) is given by 
j6 x(0 = A { x(t) + J* j5 _u ( t) + IS w(t) + F z(t) 
6 x(0) = x(0) 


(139) 


and 


§_ X.(t) = £ §_ x(t) + D 5. u(t) 


(140) 


49 



We assume that j5 _u(t) minimizes the quadratic corrective cost functional (see 
Eq. (94)) ... 


J(5 u) = 


lim E 
T-*» 




2 _'(t)Q _y(t) + 6 u’ (t)Q 6 u(t) 
-7 -ni 


+ A (t) Q*5 u(t) ]dt 


(141) 


TO j’XtriV 


Hence, the human "tracks" the nominal trajectory Xq> minimizing the deviations 
6_ 3 c(t) subject to his inherent limitations. It is assuinell that-' ' ' V ( t ) is 
"known" to the pilot so that ^(t) can be obtained from the actual observa- 
tions 2 .(t) by Eq. (137). 


The equations for determining the optimal u(t) , based on the observed 
deviations (140) , are the same as the human operator model equations presented 
in Section 3.1. 4. It is only- necessary to make the replacements 

, x(t) 6_ x(t) , £(t) + A Z(t) > «(t) + i “(t) . .. 

in Eqs. (83) - (108). The result is that the optimal regulating control component 
u_(t) is generated by 

, T^6 tt(t) + 6^ u(t) = -L _6 x(t) + v^(t) (142) 

where j5x.(t) is the best estimate of the deviations The estimate is 

generated by a Kalman filter-Predictor combination according to equations 
(102) - (108) but with X(t) replaced by _^(t) - [5. x(t) , _u(t) ] , efc. 

Closed-form expressions for the statistics of the variations ^5 x(t) , u(t) , 
j5 ^(t) are obtained from Eqs. (114) - (123) with the obvious replacements. 

Note that the variations 5_ x will contain a mean component 6_* x due to the 


The assumption T-*» is most valid when the terminal time T^ is larger than 
closed loop system time constants. T-*» conveniently simplifies the model 
structure since various feedback gains become constants. 
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"deterministic" inputs z^(t), plus a random component ^5 jc - jS x due to wind 
gusts and pilot- induced randomness. 

The total state x(0 at time t is, therefore, 

x(t) = x (t) + 5. x(t) . 


The mean, or ensemble average state trajectory is 

E{x(t) } = x(t) = 3^(0 + A x(t) (143) 

and the covariance of x(t) is 

cov [jc(t)] = cov [6^x(t)] = E{ (j$ x - x) (A * ~ A. *) ' } (144) 

By defining u^t) and the resulting x^(t) and ^(t) in the above 

manner, we have tacitly assumed that the pilot is well-trained, i.e., in an 

"ideal" situation he can control the aircraft to meet the terminal conditions. 

The decision to omit deterministic inputs z^(t) from the nominal path definition 

of x is because otherwise the control u (t) at time t would be a function 

~ ° . ~ ° [71 

of the entire future input z(o), t < O < T,. J It is unrealistic to assume 

— — — g 

that a pilot has precise a priori knowledge of z±(o ) , t £ O _< T^ especially 
since this input is a function of aircraft future attitude. By assuming that 
pilot control is dependent on only instantaneous (and past) z^t), along with 
any a priori estimates of z_(t) , the "deterministic" inputs are included 
naturally in A ■*(£)• The pilot model then estimates _z(t) continuously, along 
with the other system state variables, and uses the present estimate of z(t) 
in generating the present control. 

3.4 APPLICATION TO THE STOL LANDING PROBLEM 


In order to use the model to predict pilot performance in a landing task 
it is necessary to prespecify: 
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a. the vehicle dynamics and input disturbances, 

b. the primary display instruments, 

c. the terminal conditions, 

d. numerical values for the human response parameters, 

e. the cost functional weightings, and 

f. initial conditions at the flare window. 

For the three different STOL configurations to be studied, the vehicle dynamics 
have been developed in Section 2.4 and are presented in the state space format 

x(t) = A x(t) + 13 u(t) + F z^t) (145) 

to be compatible with the human operator model. Numerical values for A, jJ are 
given in Appendix B. Below, we discuss the display format, the a priori selec- 
tion of human response parameters and the cost functional weightings. 

3.4.1 Fully Augmented Aircraft 

With the autospeed and pitch command-and-hold systems engaged, the pilot's 
task is to control altitude using the throttles. The primary displays used are 
thus the altimeter and rate-of-sink indicator. Since we assume that a pilot 
obtains directly both the position and rate of change of displayed quantities, 
the information base is 

^(t) = col [h CG (t), h (t) , h(t) ] (146) 

where h (t) is obtained from the altimeter and h(t) and h(t) from the rate- 
* 

of-sink meter. With the aid of Eq. (81) displayed outputs may be written as, 


* • 

h(t) is also obtained (redundant information) from the altimeter — but not to 
the degree of accuracy as from the fi instrument directly. 
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t > 0 


(147) 


2 .(t) = C x(t) +2o(t), _ 


where, for the augmented aircraft, x = col [C g , <*' , 6h', 6T] and 


C = 


0 

0 

7.9 


0 

-126 

70 


126 

0 

0 


0 

0 

16.9 


(148) 


The additional term ^ (t) = col[53. - 13. 2t, -13.2, 0] is included since the 
mathematical state x represents deviations from an equilibrium 6° glideslope. 

Figure 4 shows the flare geometry. The flare maneuver begins at an altitude 
h = 41 ft. (h„_ =53'). At this altitude, on a 6° glideslope, the plane is 
at R = -140' from the runway threshold. The desired touchdown point is taken 
as the center of the touchdown zone. Thus, R^, =475. Since the airspeed 

is V q = 126 ft/sec, the nominal touchdown time is 

T d -M ' 4 - 9 sec - 

In the mathematical model, Eq. (145), 126x^ = h is the altitude deviation from 
the equilibrium 6° glideslope. In order that h(T^) =0, we require 

<$h(T ) = 225 tan 6° = 23.53 (149) 

a 

Sink rate at touchdown should be no greater than -3 ft /sec. We, therefore, 

. 

choose a nominal h(T,) = -2 ft/sec. Since h = -13.2 + 6h, the boundary 

. d 

condition on 6h(T,) is 

d 

<$h(Tj) = 11.2 (150) 

cl 


The boundary conditions (149) - (150) may be put in the form 
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FIGURE 4 TOUCHDOWN GEOMETRY 
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(151) 


H x(T) + c = 0 
a — 

where 



'o 

126 

0 

0' 


'-23.53' 

H = 





> £. ° 



.0 

0 

-126 

0 . 


.-H.2 . 


Nominal values for the human response limitations are chosen on the basis of 
past experience with the optimal control model. As discussed in Section 3.1, 
these values are 

T - . 2 sec 

p = .003 tt = -25dB (152) 

u 

P yl = P y2 = P y3 = * 0lTr /- 5 ° - 17dB 

where the observation noise/signal ratios have been adjusted to -17dB in assuming 
equal allocation of attention to both instruments. Any effects of pilot monitor- 
ing other instruments (e.g. , pitch attitude, engine RPM) that are not needed 
explicitly for control are neglected. Thresholds, a^, for the instruments are 
dependent upon display gain and scale markings, and relate to the accuracy 
with which one can read a given variable. For the standard type set of instru- 
ments used in this study, a priori values 

2 

a^ = 5 ft, = .8 ft/sec, a^ = .4 ft/sec (153) 

were chosen rather crudely. a^ and a ^ correspond to one-half the minimal scale 
division on the altimeter and rate-of-sink instruments, respectively. By assum- 
ing that a signal must move through its threshold in 2 sec in order for motion 
to be discerned, we obtain = .4 ft/sec . Fortunately, model predictions 
were not found to be highly sensitive to changes (by a factor of 2) in the 
threshold values. Thus, for the STOL landing problem, the numbers a. can be 
assumed to be (within a factor of 2) representative of the standard h, h 
instruments. 
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The cost functional weightings in JOS u) are determined a prior as dis- 
cussed in Section 3.1.3. These weightings relate to the task of optimal regula- 
tion about the nominal flare path. Representative values can be determined in 
part from analyzing the approach task, which is totally a regulation task. We 
first choose a weighting on altitude deviations 6h. During the approach, it 
is assumed that the pilot desires to maintain glideslope error within +.3°. We 
assume, therefore, 

Sy - .3° = .0052 rad 
max 

Since Sh = RSy, we have at the flare initiation point, 

6h = 390*(. 0052) = 2.05' . 
max 


This value of 6h ma y is used for the entire flare path. The weighting 
is, therefore, by Eq. (95), 


q 6h " 


1 

(2.05) 2 


.25 


(154) 


Since the throttle input is limited in magnitude, ST c is weighted in the 
cost functional. On the approach 5T c = .98 which corresponds to 13,200 lbs of 

thrust; 6T increases nominally to 1.4 during flare. Since the range of ST 

c 2 C 
is approximately 0 - 2.4, 6T , max - 1.0. Therefore, the weighting on <$T is 

c c 

chosen 



1.0 


(155) 


The final weighting to be chosen is the control rate weighting. As discussed 

in Section 3.1.3, q^ is chosen to yield a - .1 sec. For the augmented 

c 

aircraft, the required weighting was found to be 

q 6T = -022 (156) 

c 
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This corresponds to a maximum rate of throttle motion of -7 units of thrust/sec, 
or movement through the entire throttle range 0 - 2.4 in approximately .3 sec. 

Finally, the initial conditions on the state covariances appropriate to 
flare initiation must be specified. These were obtained from the model by 
"flying" the aircraft down the approach path and stopping at an altitude of 
h^^ =53'. The predicted variances of the system state variables at h = 53' 
were then used as initial conditions for the flare segment of the landing task. 
The values used were 

°Sa' = -005, cr^, = *006, cr^ = .06, cr^ = *1 (157) 

c 

These values were subsequently compared with NASA simulator data taken at the 
flare initiation window. Generally, model and data results agreed within a 
factor of 1.5. Slight errors in these quantitites are relatively unimportant 
as regards touchdown, since the effects of the initial model covariances damp 
out after approximately 1 sec. Initial conditions on the mean state x.(0) are 
zero, since the plane is assumed to be on the glideslope (on the average) and 
the C ground effect has not as yet built up. 

L* 

3.4.2 Aircraft with AutOspeed Only 

With the autospeed system, the aircraft can approach the flare initiation 
point in the 2° nose-up attitude desired for touchdown. As discussed in 
Section 2.3.3, the flaps fif^ must be backed-off -13.1°, and the equilibrium 
tail position is fit = -3.5°. The pilot is required to regulate pitch devia- 
tions from 2° through controlling the tail input. Pitch deviations occur 
through cross terms in the throttle input, and directly from the ground 

effect. 


This is necessary to avoid transients in the model response over the first 
second of the flare. These might arise if the model were started "cold". 
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The equations of motion for the STOL dynamics with autospeed are given in 
Section 2.4.5. For this task, the pilot's major displays are the altimeter and 
rate-of-sink meter as well as the pitch attitude indicator on the eight-ball. 

The pilot's information base is, therefore, 

£(t) = col [h CG (t), h(t) , h(t) , 0(t), q (t) ] (158) 

or ... <•; . 

^_(t) = <3 x(t) + D u(t) + ^(t) t _> 0 , . . (159) 

where in this case x = col [C^, C^, a', q, 0, 6h' , ST] , u = col [ST c> St] and 
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The equation for h(t) is obtained by substituting for 0 - a' = Sh'. 

The terminal conditions on Sh and Sh for this problem are given, as 
before, by Eqs. (149) - (150). In addition, terminal conditions are placed on 
S0 (deviation from 2 ° equilibrium) and q, i.e., 

0(T d ) = 0., q(T d ) = 0. (160) 

The four terminal conditions are put in the general form (151) with 
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Nominal values for the human response parameters are 


T = .2 sec 


p = p = . 003 TT = -25dB 
u, u„ 


TT = ~ “ 15dB, i=l , 5 


(161) 


The observation noise/signal ratios have been adjusted to -15dB assuming equal 
allocation of attention to all three display instruments. The thresholds a^ 
for the altimeter and rate-of-sink indicator are given by Eq. (153). For the 
pitch attitude indicator we assume, 

a^ = 1.0 deg , a,. = .5 deg/sec. (162) 

The cost functional weightings for the autospeed case must reflect the fact 
that the pilot is required to regulate both altitude deviations and pitch devia- 
tions from the nominal flare path. The cost functional weightings on 5h and 

<ST are chosen as in the augmented case. Thus, 

c i > 

q 6h = . 25 q 6T =1. 0 ' (163) 

c 


In the past modeling efforts using the optimal control model, pitch regula- 

[3-4] 

tion was treated by including a cost functional weighting on q. This 

assures that pitch error corrections are made smoothly, and corresponds with 
pilots' subjective control behavior. In the present study, we assume 
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6q =1 deg/sec 
max 

Thus, the cost functional weighting on pitch rate deviations is 

% = (164 > 

For the STOL aircraft under consideration, it was found that weighting 
q(t) alone resulted in large pitch angle deviations during both approach and 
flare. This is due to the poor pitch damping of the aircrafts 1 It is assumed 
necessary for the pilot to control pitch angle directly; we assume a maximum 
tolerance 


60 = 2 ° 
max 

so that the aircraft would land with 0^0. Hence, 

■•«9 '(I ) 2 "; 25 1 (16 

The control weighting on tail deflection 6t is determined from the tail 
position limits. The tail is free to move through a total range -10 to +10 
degrees. Since the trim position is -3.5°, the maximum allowable deviation is 
6.5° - .11 rad. Thus, 


q 6t = 



(166) 


The control rate weighting on commanded throttle u^ = 6T c was chosen to 
yield a resultant T = .1 sec. The required value was found to be (as in 
the augmented case) 



c 


(167) 
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The control rate weighting on tail position = 6t was chosen to yeild a 

The required control rate weighting q • was found 


resultant 
to be 


x . = . 1 sec. 
n2 


6t 


- 16. 


(168) 


which corresponds to a maximum tail rate of approximately .25 rad/sec - 15°/sec. 

Initial standard deviations for the model states were determined from 
prediction at the flare window. For the autospeed case, it was found 


6a' = .008, 6q = .003, 68 = .006, 6h’ = .003 

6t = .14, 6t = .23, 6t = .004 
c 


(169) 


One simplifying assumption is made in modeling the pitch ground effect in the 
autospeed case. The ground effect begins at an altitude of h^ = 60'. 

We wish to start the flare model at h = 41'. Prior to flare initiation, 

LCj 

the pilot's major role is to regulate pitch deviations, including those that 
arise from C between h = 60 - 41'. We, therefore, assume that pitch 

M Lvj 

deviations 60 = 6q = 0 at h = 41', so that it is only necessary to consider 
the additional ground effect that is introduced for h^ _< 41. This 

enables us to begin the model propagation at h = 41' (t=0) with 6x(0) = 0. 

• • 

The equation for generating C^(h(t)) remains as in Eq. (76), 

but with 0^(41') = 0. The additional pilot tail input at h = 41' is computed 
from the average value of between h = 60' and 41' which is - .05. 

From the q equation, the tail deflection necessary to balance this average 
moment is -.5°. Thus, since the equilibrium tail position is -3.5°, we have at 
flare initiation 6t - -4°. 

3.4.3 Unaugmented Aircraft 

The basic aircraft approaches flare initiation in a -4.1° nose-down con- 
figuration. The pilot must rotate the plane to achieve a 2° nose-up configura- 
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tion at touchdown. Since the aircraft is pitched down, a higher equilibrium 
thrust (17,500 lb) is needed to generate the required lift to maintain equilib- 
rium. 


The equations of motion for the basic STOL dynamics are given in Section 
2.4.5. For this task, as for the autospeed case, the pilot's main displays 
are the altimeter and rate-of-sink meter, and the pitch attitude indicator. 

The display information base is, therefore, given by Eq. (158). This may 
be written 

2 _(t) = £ x(t) + D _u(t) + ^(t) (170) 


with x(t) 

= col 

[ <V 

C M’ U ' 

, a’, q 

9, 6h', 

St], 

and 
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£ (t) = [53 -13. 2t, -13.2°, 0, -4.1, 0]. 
o 

The equation for h is obtained from differentiating <5h' = -.105u' - a' + 0 
and subsituting for u' , a' , q. 

The terminal conditions on <5h and 6h are given by Eqs. (149) - (150). 
The terminal conditions on 9 (deviations from -4.1°) and q are 

9(T d ) = +6.1, q(T d ) = 0. (171) 
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In the general form (151) we have 


r. 
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For the basic aircraft, we do not assume that the pilot controls speed 
directly (indeed speed variations on the flare are small). Thus, we have not 
included u' as a displayed variable, nor is it necessary to include a weight- 
ing on u' in the cost functional. The nominal human response parameters for 
this case are, therefore, the same as Eq. (161). The thresholds a^ are the 
same as those before, namely 

2 

a^ = 5 ft, a 2 = .8ft/sec, a^ = .4 ft/sec , = 1.0 deg, a^ = .5 deg 

The cost functional weightings on 5h, 6q and 60 are the same as those 

chosen for the augmented case. The cost functional weighting on 6T is 

. c 

somewhat higher than before since the available throttle range from nominal 

■k 

is smaller than in the autospeed and augmented case. Thus, 


q 


6T 

c 


1.5 


(172) 


Also, since the equilibrium tail position is -1°, ( ^ t max = 9°. Hence, 



(173) 


The control rate weightings on throttle and tail are adjusted to yield 


T = 
n. 
x 


. 1 . 


* 


nominal throttle setting is 17,500 lbs as opposed to 13,200 lbs. 
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The required values are 



.033, = 16. 


(174) 


The initial standard deviations for 6. x were determined from model predictions 
at the flare window. The values found were 


6u' = .02, 6a' = .01, <5q = .003, 66 = .006 


6h' = .01, 5T = .14, 6T = .23, St = .004 

c 


(175) 


As in the autospeed case, the pilot was assumed to have compensated for the 
ground effect over the range h = 60 - 41' so that the model could conveniently 
be initialized at t=0 at h = 41'. This introduces an additional tail deflec- 
tion at flare initiation of .5°. Thus, initial tail deflection is -1.5°. 


3.4.4 Model Initial Estimates 


It is necessary, when applying the optimal control model to study problems 
involving "deterministic" inputs, to specify the human's a priori estimate 
x(0) of the system state, x, as in Eq. (107). For the landing task being 
studied, the aircraft is nominally on the glideslope at flare initiation. Hence 
vehicle states which represent deviations from the equilibrium 6° glide-path 
are zero at t=0, and it is reasonable to assume that the pilot's estimate of 
these states is also zero. Therefore, it remains necessary to specify the 
pilot's a priori estimate of the ground effect terms x^ = 6^, = G^. 

A A 

We assume values for the initial a priori estimates C^, equal to 20% 

of the maximum values for C T and C„. 

L M 


Thus, 


C M (h LG * 60,) " -° 6 


C L (h LG * «’> * -° 9 
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As discussed in Section 3.4.2, the initial estimate on C w enables us to intro- 

M 

duce a simplification in the model. Since -C^ increases from 0 to .1 
between h = 60' and h = 41' , it can be assumed that the pilot, because of 
his initial estimate, has corrected for the pitching moment during this 

time interval. Thus, we assume 60 = 0. and 6q =* 0. at h = 41', and it 
becomes necessary to consider only the increase C = .1 for h _< 41'. With 

A 

an initial estimate = -.06, the pilot's estimate of by the time 

h = 41' will be approximately -.1 + .02. Thus, it can be assumed that the 

pilot's estimate of the increase C + .1 at h = 41' is zero. 

M 

As a result, the flare portion of the flight including the pilot reaction 
to the ground effect can be modeled as beginning at h - 41', where we set 
t=0 , with 


= -.09; $ 2 = 0; x ± = 0, i=3,...,n 

where + . 1. 

3.5 TOUCHDOWN PREDICTIONS 


The independent variable in the pilot model is time. At any time, t, the 

model generates predictions of the mean and standard deviation of all system 

variables according to the covariance propagation equations derived in 

Section 3.2. In addition, as discussed in Section 3.2, these quantities are 

used to generate the state and output probability densities as functions of 

time. Of prime interest in a landing task are the probability densities and 

mean values of key system variables (especially sink-rate) at touchdown . We 

have model predictions available as functions of time . The transformation from 

the time predictions to the touchdown predictions is complicated because the 

touchdown time, T, , is a random variable (E{T,} =: 4.9 sec). Below it is 

d a 

shown how we can generate touchdown predictions, including the probability 
densities of touchdown time (which translates into range dispersions) and sink- 
rate. 
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For ease of exposition, let z^(t) = 'h G (t)' and let z 2 (t) be any model 
output whose mean and/or probability density at touchdown is sought. . Let 


nu(t) = E{z^(t) } 


i=l,2 


(176) 


cr (t) = E{ (z (t) - m )(z.(t) - m.)} i,j=l,2 

1 11 


• fut-vn n-orittutoo -.-u- 


be the mean and covariance at time t of the Gaussian random variables z^ and 

z„. The quantities m. and o. . are available directly from the results of 
2 i ij , 


Section 3.2. 


'f : : s 


The joint probability density p(z^(t), z 2 (t)) is written as 


pU 1 ,z 2 ) = 


- p exp - 

27r[a ll d 22 (l - p )* 


Z ! - 2Z 1 Z 2 + Z 2 


2(1 - p A > 


(177) 


J 


where 


Z. = 


h - m i (t) 


V° ll (t) ’ 2 V 0 22 (t) 


z 2 - m 2 (t) 


(178) 


and 




(179) 


is the correlation coefficient between z^(t) and z 2 (t). The univariate 
probability density of z^(t) is 


p(z (t)) = — — 

x y 2ou 


exp 


4 


(180) 


11 
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Let be the touchdown time and let p(T^), to be determined, be the prob- 

ability density of T^. The mean of z 2 at touchdown is given by 


■/' 


E(z 2 (T d )} = J E{z 2 (t)|t = T d > p(T d ) dT d 
0 


where the condition mean is defined as 


(181) 


2 1 d 


■/ 


E ^ z o T.} = / ZoP(Zo T.) dz 


2 1 d 


(182) 


In the process of integrating Eq. (181), we do not wish to include those trajec- 
tories that overflare. The model that has been developed is not appropriate for 
predicting pilot response to an overshoot of the landing zone. Thus, in comput- 
ing E{z 2 (Tj)} from the model predictions only those flights that land within a 
prescribed time T - T^ may be considered. Flights that have not landed by t= 
are assumed to have overflared. In our analysis, T is taken to be the first 
time t _> T^ = 4.9 sec at which E(h(t)} = 0. The result of considering T^ < T 
is to modify Eq. (181), i.e. , 


E{z 2 (T d ) |T d < T} 


f 


E(z ? (t) jt = T d )p(T d )dT d 


/ 


P< T d> dT d 


(183) 


It thus remains to compute E { z 2 J T d } and p( T d )* By definition, the condition- 
al density p ( z £ | ) is also given as 

p(z 2 (t)|t = T d ) = p(z 2 (t) |z^(t) = 0) (184) 
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since z^(t) = 0 implies t = T^. The Gaussian density p (z^ (t) | (t) is 
obtained by dividing Eq. (177) by Eq. (180), yielding 


P(z 2 l z i> = 


[ 2 tto 22 (1 - p 2 )]' 2 


exp 


(z 2 - P z x ) 2 

2(1 - p 2 ) 


Evaluating the conditional mean (182) with z^ = 0 gives 


E{z 2 (t) |t = T d > = - p 


The density of T^ is obtained by noting that 


°22 a i2 

— mi + m 2 = - — u»i + m 2 


11 . 


11 


. (185) 


Pr (landing between time t and t + dt) 
= p(T d )dt 
= -p(h(t) = 0)dh 


Thus, the required density is 

p(T d ) = - p (h=0) ~ = - p (h=0) • E{h(t) } 


(187) 


where 


p (h=0) = 1 e m i/2a (188) 

/2ua^ 

and where E(h(t)} = the mean sink rate at time t is obtained from the model 
predictions of Eqs. (114) - (118). Therefore, the mean of z 2 at touchdown 
may be computed from Eq. (183) upon substituting Eqs. (185) and (188). Note 


68 



that it is possible to carry through the above process to compute any moment 

E{zo(TJ} n > 2, in addition to the mean z.(T ). 

L d — 2d 

The probability density of z~(T,) in general can be computed by first 

L d 

finding p(T |z ) and then using Bayes rule, or by numerically obtaining a 
histogram. Either method is cumbersome. Analytic results are obtainable when 
E{z 2 (t)} is a monotone function of t, i.e., El^Ct)} has the same sign 
for all t > 0. Sink rate z_(t) satisfies this condition. Consider a given 
time t. From Eq. (185) we can evaluate E(z 2 (t)|t = T^} = z; we wish to 
'determine p(z(T^)) evaluated at z(T d ) = z. The probability that z(T d ) lies 
between z and z + dz is equal to the probability that touchdown occurs 
between time t and t + dt since z^(t) is monotone in t. Thus, 

p(z(T d )) = p(T d ) 

z(T d ) = 2 

Eq. (189) will be used to compute the probability density of h at touchdown. 



69 



IV. MODEL RESULTS 


In this section, model covariance propagation results are compared with 
experimental data obtained on the NASA-Langley Research Center Real Time Dynamic 
Simulator. Three cases are studied; namely, the fully augmented STOL, the basic 
aircraft with autospeed and the unaugmented aircraft. Thus, we will be able to 
examine the effects on pilot response of the two augmentation schemes. 

4.1 DESCRIPTION OF EXPERIMENTS ~ ; 1 


The experiments were conducted on the NASA Langley "Research Center Real 

• • «_• - 1 j ' 

Time Dynamic Simulator. The simulation was fixed-base, and only the longitud- 
inal flight model was simulated on the CDC 6600 digital computer with table 
look-ups used for the various aircraft coefficients and throttle/thrust charac- 
teristics. 

The aircraft was positioned initially R = -5000 ft from the runway thresh- 
old. The pilot had to acquire the glides3,ope, and then trim the aircraft to 
maintain equilibrium flight down the 6° glideslope. Flare initiation was at 
h = 53', at which point a flare initiation light came on. A "get ready" light 

CLi 

came on at approximately 5 sec prior to flare initiation. The flare geometry 
is shown in Figure 4. The pilot was instructed to land the aircraft, and no 
go-arounds were permitted. Only the basic set of instruments were available 
for landing the aircraft: pitch attitude indicator, altimeter, sink rate, and 
airspeed, dials. No external visual cues were provided. 

4.1.1 Data Analysis 

Ten flights of each of the three conditions .... 

• Fully Augmented Aircraft 

• Autospeed System Only 

. • Basic Aircraft , 
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were made in sequence. Flare (and approach) data consisted of the sampled 
time histories of all pertinent system variables to touchdown. The sampling 
time was .25 sec during flare. Data from only one pilot was provided for 
validating the model developed in the present effort. 

Ensemble statistics were computed for various flight quantities during 
flare. For a given variable, x, the sample mean and standard deviation as 
functions of time were computed according to 


X< V ' I z W 

k=l 


(190) 


a 

X 


<'i> 


1 

N 


i 

N 

£ 

k=l 


x k (t i ) ' 


(191) 


where x^(t) is the k 1 "* 1 sample path and where N is the total number of sample 

paths (N = 10). For convenience, we define t = 0 to be the time at which the 

aircraft is at h =53'. 

CG 

For the augmented STOL, the ensemble mean and RMS were computed for sink- 

rate (h) , throttle (T ) , total thrust (T) . For' the autospeed and unaugmented 

c 

cases, additional statistics were also computed for pitch angle (0) and tail 

deflection (St). In the- following section, these human-generated statistics 

are compared with the model predictions. However, it should be borne in mind, 

when viewing the comparisons, that because of the relatively low number of sample 

. 

paths averaged, the variance in x(t.) and O (t) is high. Therefore, these 

X X 

experimental quantities can only be viewed as approximations to the true en- 
semble statistics. 


— /S 

The variance in the sample mean and RMS decreases as N . Thus, with N = 10, 
any one sample path can have a pronounced effect on the resultant ensemble 
statistics. 
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4.2 MODEL-DATA COMPARISONS 


In this section, model covariance propagation results, obtained using the 
expressions derived in Section 3.2, are compared with ensemble statistics com- 
puted from the pilot-simulator data. In addition, probability densities of 
touchdown time and sink-rate at touchdown (Eqs. 187 - 189) are presented. The 
model predictions have been generated using the model parameters as chosen 
a priori in Section 3.4. Thus, the model results are a priori predictions, and 
not the consequences of model-data matching. After presenting'-'the results for 
the three different STOL configurations, model attributes and deficiencies will 
be discussed indicating possible changes in the a priori parameter values and 
the model structure. ' 

4.2.1 Fully Augmented STOL 

With both the autospeed and pitch hold systems engaged, it might be expected 
that the pilot is capable of controlling the aircraft with minimal flight path 
dispersions. Figure 5 shows h (t) for each of the ten data trails, i.e., a 
scatter diagram of vs - fc * The dashed lines enclose the model's pre- 

dicted 95% confidence limits on altitude dispersion. 


Ah(t) 



= [h(t) - 20 h (t), h(t) + 20 h (t)] 


(192) 


where h(t) and 0, (t) are the model predictions of mean altitude and variance 
** h 


at time 


For t > 4.9 sec, which is the nominal touchdown time, model 


predictions are obtained by setting u (t) = u (4.9) = constant, since 

nom nom 

u (4.9) = 0. 
nom 


Indeed, Systems Control, Inc. did not receive the NASA data until after model 
predictions were generated. 

Since h(t) is assumed to be a Gaussian random variable, the 20 band about the 
mean is the 95% predicted confidence limit. 
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As can be seen in Figure 5, model and data dispersions are in general agree- 
ment, with the pilot showing a somewhat greater tendency to overflare. Three 
flights have not yet terminated by t = 6 sec, and one trajectory is on the 
outer limits of any distribution. The probability density p(T^) of touchdown 
time (Eq. (187)) for T^ _< 6.5 is superimposed on the dispersion diagram. This 
density is non-Gaussian. p(T^) increases rapidly for t > 4 sec and has a 
well defined peak between 4.25 < T^ < 5.5, indicating a high likelihood of 
touchdown in this time interval. As t increases, p(T^) decreases more slowly, 
indicating a somewhat greater tendency to land long as opposed to landing short. 
Because of the low number of sample runs, a histogram of T^ was not computed 
for comparison with the predicted p(T^). 

A 

The probability of landing prior to t = T sec is given by 


Pr (T J < T> 
d — 



A 


T 

p( T d )dT d 


0 


A 


P(T) 


(193) 


The value of T such that P(T) =0.5 is defined as the median touchdown time, 

* 

and must coincide with the touchdown time ; of the mean flight path h(t). For 

this case the median touchdown time is 5. 15 sec. We note that this is not equal 

to the average touchdown time due to the asymmetry of p(T^). Evaluating the 

integral (193) for T = 5.9 and T = 6.5 yields P(5.9) = 0.69 and P(6.5)=0.73, 

respectively. This indicates that approximately 30% of all flights are predicted 

to land more than one second later than the nominal T^ = 4.9, and approximately 

25% land more than 1.6 seconds late. ' Translation into range at touchdown 

(R = V T,) implies an approximate 0.25 probability of overflaring the touch- 
td o d 

down zone. It should be emphasized again that model predictions are not valid 
much beyond t = 6.0 sec. Once a pilot realizes that he is about to overflare 
his (cost functional and) control strategy may change drastically, e.g., diving 
for the deck, cutting all power, etc. Thus, although the model predicts that 
30% of all flights land at T^ > 6.0 sec, the existing model cannot predict the 
subsequent pilot control inputs. 


Since the density of h(t) is Gaussian, it is symmetric about its mean. 
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Figures 6 and 7 compare model and data sink rate, throttle input and result- 
ant thrust. Model predictions consist of the mean (solid line) and the +1 o 
band about the mean (dashed line). The pilot data consists of the ensemble 
mean (Eq. (190)) and the +lcr ensemble deviation (Eq. (191)) at .25 sec intervals. 
The model and pilot data are in reasonable agreement. The pilot's mean response 
generally lies within the 1-sigma band of the model predictions. The model and 
data standard deviations agree generally within a factor of two. Model and pilot 
sink-rates agree within 0.5 ft/s ec over the entire flare. Note that the average 
pilot sink-rate at t = 0 is -13.6, which is slightly greater than the assumed 
nominal of r -13.2. • 

Both model and pilot throttle inputs increase rapidly during the first 

.5 sec. This is primarily in response to the initial estimate of the C ground 

Li 

effect. The pilot's commanded T continues to increase for .5 < t < 1 

c 

whereas the model does not show this trend. Possible reasons for this slight 

discrepancy include the slightly greater pilot initial sink-rate as noted above, 

an initial pilot estimate of C greater than the -.09 assumed in the model, 

1 • * 

(or an a priori pilot estimate of both C and an initial rate of change C ) , 

1* L 

or slight differences between the actual throttle/thrust characteristics and 
the simplified first-order lag assumed in Eq. (66). On the average, the net 
throttle input of the pilot is slightly greater than the model's, yielding a 
pilot sink rate approximately .5 ft/sec lower than the model as touchdown 
approaches. 

The basic assumption in modeling the terminal time aspects of the landing 
task was that in the absence of all external disturbances, the pilot would 
generate a minimum rate open-loop control input that met the terminal state 
conditions. In the presence of external disturbances, the pilot would then try 
to regulate about the so-generated nominal flight path. The pilot data provides 
an opportunity to validate this hypothesis. Figure 8 shows the sink-rate and 
throttle for the nominal flight path (solid curve), mean model predictions that 


Section 4.3.1 discusses the effects on model predictions of other characteriza- 
tion of C L , specifically allowing the pilot model to continually estimate both 
Cl and its rate of increase. i 
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FIGURE 6 SINK RATE, AUGMENTED STOL 
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FIGURE 7 THROTTLE/THRUST INPUTS, AUGMENTED STOL 
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FIGURE 8 REGULATION ABOUT NOMINAL FLARE PATH, AUGMENTED STOL 
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include all external disturbances (dashed curve) and the pilot data ensemble 
average. There is an excellent agreement between the mean pilot sink rate and 
the nominal sink rate. This indicates that our assumption of pilot regulation 
about a minimal control rate nominal flight path is most reasonable. Indeed, 
the pilot seems to do a better job than the model in regulating disturbances 
about the nominal altitude path. This is probably due to the pilot's having a 
better internal characterization of C than was assumed in the model. The 

Lt 

difference between the nominal T and the model/pilot T in Figure 8b is 

c c 

the additional control^ input T applied to correct for C . The pilot appears 
to have corrected more precisely, and earlier, than did the model. 

4.2.2 STOL with Autospeed Only 

With the autospeed system engaged, the pilot can still maintain a +2° pitch 
attitude on the approach. During flare, it is only necessary for the pilot to 
regulate the deviations in 0-2°. Consequently, it might be expected that the 
altitude control (i.e., flare) behavior of the pilot would not differ substan- 
tially from that observed in the augmented case. Figure 9 shows the altitude 
dispersions from of the sample flare paths, and the 95% model confidence band. 
Comparing the pilot data in Figure 9 with that of Figure 5 (discounting the one 
or two "strange" runs) shows a noted similarity in altitude dispersions in the 
vicinity of touchdown. The tighter clustering of the scatter points in Figure 9 

indicates that the pilot's control, strategy has probably become less variable 

* 

by the time the autospeed experiments were performed. Therefore, external 
effects of learning may be biasing the experimental results since the pilot's 
control strategy is continuing to stabilize on a run-to-run basis. , 

Comparing the model flight path statistics of Figure 9 with those of 
Figure 5, we see an almost identical variance but a mean flight path that lands 


The augmented STOL runs were the first series of 10 experiments, and no warm-up 
runs were made. The autospeed case was the third set of 10 experiments actually 
performed. An intermediate augmentation scheme involving pitch damping was also 
simulated. However, this case is not discussed herein. Thus, a possible ex- 
planation of the decreased variability is that by the time the autospeed-only 
runs were made, the pilot had "warmed-up". Obviously, this is a short-term 
adjustment phenomenon only. 
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FIGURE 9 FLIGHT PATH DISPERSIONS, AUTOSPEED STOL 


sooner at t = 4.7 sec. The scatter points of Figure 9 also indicate that our 
prediction of altitude dispersion is reasonable, but the model is slightly low 
in its mean altitude and touchdown time predictions. The probability density 
p(T^) is also shown in Figure 9. The results are much the same as in Figure 5, 
but with a .5 sec shift towards earlier touchdown times. Indeed, integrating 
Eq. (193) with T = 6.5 indicates that only 20% of all flights might be expected 
to overflare compared with 25% in the augmented case. Unfortunately, the trend 
towards lower touchdown times does not appear in the pilot autospeed data. 

Comparisons of model vs. data statistics of sink-rate, throttle input and 

thrust are shown in Figures 10, 11a and 11b, respectively. In all cases, the 

ensemble standard deviations are much smaller than their counterparts in 

Figures 5-6. This is further evidence of the reduced variability in the pilot 

response caused in part by underlying adjustment phenomena inherent in the 

* 

experimental procedure. The model predictions of the standard deviations a^, 
a Tc and o^, are noticeably greater than the pilot data. Several reasons for 
these discrepancies are: 

a. The pilot may be devoting more attention to the sink-rate indicator 
than we have assumed. More precise estimates of h will result in 
less throttle variability and more accurate h control. 


Cost functional weightings, other than the ones chosen, may be more 
representative of pilot behavior. For example, weighting = <$h 
or y^ = <$h in the cost functional would have a very pronounced 
effect on reducing both throttle and sink-rate variability. 


c. A lower value of the motor-noise ratio pu may be appropriate. This 
will result directly in less control input variability, and, hence, 
less closed-loop randomness. 

Indeed, with "equal" training in both augmented and autospeed landing tasks, 
one would expect the pilot to show a greater variability in the latter since 
it is a somewhat more difficult overall ^:ask. It is implicitly assumed that 
the pitch-hold augmentation is working properly (i.e. , q - 0) , and that adverse 
interaction between pitch-hold augmentation and pilot control is minimal during 
flare, allowing, in effect, pitch-loop decoupling. 
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FIGURE 10 SINK RATE, AUTOSPEED STOL 
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FIGURE 11 THROTTLE/THRUST INPUTS, AUTOSPEED STOL 


d. A larger value of T n may be appropriate for the specific pilot. This 
would lower the bandwidth, and, consequently, the power, of the injected 
motor noise (see Eq. (97)). 


Items (c) and (d) suggest that a reexamination of the model's treatment of 
motor noise is in order. In fact, other applications of the human modeling tech- 
niques to study time-varying problems have arrived at similar conclusions. ^ 

Comparison of the model mean with the ensemble average for h, shows a 
well-defined difference between model and pilot response. For 1.5 < t < 3 sec, 
the model's sink rate is considerably (.3 - .8 ft/sec) higher than the pilot's. 
This difference arises from the higher level of pilot thrust (yielding more lift) 
between .75 < t < 2 sec as shown in Figure lib. For t > 2 sec, the model 
must correct the compounded effects of the high sink-rate and more thrust is 

"k 

applied to put the aircraft back onto its nominal flare path. The net result, 
as seen, is that the model's mean flight path lands somewhat earlier than the 
nominal. The cause of the model's behavior can be traced directly to its esti- 
mate of the C ground effect. The model corrects for its initial estimate 
-L 

C = -.09, but does not anticipate a further increase in C until 1 sec 

L ii ^ ## 

later — when the effects of the C T are eventually observed in h, 6h and 6h. 

L,ge 

Correction for C is, therefore, applied late. The pilot, on the other 

L,ge 

hand, seems to be aware earlier of the ground effect. He compensates earlier 
and needs to generate less lift nearer the ground. 


The comparisons of both the average pilot response and the mean model 
response with the nominal sink rate are striking. Figure 12 shows that the 
pilot is doing an excellent job of regulating h about the open-loop model- 
generated nominal. Our basic assumption that the pilot regulates about a "least 
workload" (minimum u) nominal flight path is demonstrated convincingly for the 
second time. However, our model for how this regulation task is accomplished, 

and how. the pilot compensates for C needs modification. In Figure 12b, 

I* j 8® 


We see here the interplay between the open-loop generation of a "nominal" flare 
path, and the closed-loop requirements to minimize deviations from this flight 
trajectory. 
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FIGURE 12 REGULATION ABOUT NOMINAL FLARE PATH, AUTOSPEED STOL 


the additional throttle inputs ST needed to counterbalance C„ indicate that 

c L 

the model requires a larger initial estimate of C , or an estimate of C 

Li Li 

that increases with time, to generate a control input that would match the pilot's. 

The aspects of the pitch regulation task are shown in Figure 13. It is 
apparent that the pilot has devoted little, if any, attention to regulating 
pitch about the 2° nominal. The elevator input is essentially held constant 
during the flare. Pitch excursions are highly variable. The model predictions 
also show virtually no elevator input to regulate 0 which is predicted to 
decrease to 0.5° at touchdown. The major difference between pilot and model 
pitch response lies in the variance. The most plausible reasons why pilot vari- 
ability is larger than the model's are the following: 

a. The pilot is probably devoting most of his attention to the altimeter 
and/or rate-of-sink instruments. Therefore, the observation noise 

on 0 and q would be higher than the -15dB assumed, resulting in 
poorer estimates of pitch deviations and, consequently, less control 
of them. The fact that pilot variability of h, h are lower than 
the model, while variability of 0 is greater than model predictions 
gives further evidence to greater pilot attention to altitude control. 

b. The pilot's strategy in controlling 0 may show a run-to-run vari- 
ability. Indeed, the sensitivity of model predictions with respect to 

2 

changes in the cost functional weighting on 60 is high and can 
account for the observed variability. The nominal weighting on qg 
is .25. Varying q Q over the range .25 +.25 results in pitch angles 
at touchdown between 0.5 +1.0°, yet has little effect on other model 
predictions. The high sensitivity of 0 to small variations in pilot 
response (caused by the poor inherent pitch damping of the STOL) is 
one reason for the poor pilot rating. 


In past applications of the optimal control model to pitch regulation tasks, 
only pitch rate was included in the pilot's cost functional. The fact is that 
with qg = 0, pitch deviations are -2.5°, which indicates a serious handling 
quality deficiency. Thus, the pilot must adopt a pitch regulation strategy 
somewhat different from his "usual". 
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FIGURE 13 PITCH CONTROL RESULTS, AUTOSPEED STOL 
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4.2.3 Unaugmented STOL 


The model predictions and pilot data for the unaugmented STOL are shown in 
Figures 14 - 17. From the results, it is apparent that: (a) the aircraft is 
near impossible to fly, or (b) the pilot is not well-trained on this task. The 
fact that the model seems to control the aircraft is evidence of the latter. 

The flight path dispersions are shown in Figure 14. Only eight trajectories 
are shown, yet it is apparent that pilot variability is exceptionally high, 
with a marked tendency to overflare. The model predicts a higher variability 
in altitude dispersions for the basic STOL configuration than in the jpreceeding 
cases. The model's median touchdown time, however, is only 4.65, i.e. ; the 
model shows a tendency to land early, as in the autospeed case. The probability 
of landing at t > 6 sec is only 0.2; this seems low in view of the observed 
pilot trajectories. 

Sink-rate and throttle inputs are shown in Figures 15-16, respectively. 

Pilot variance is extremely high, about twice as great as the model predictions. 
The large variance leads to the belief that pilot strategy is highly variable 
on a run-to-run basis, so that the effects of learning have not yet stabilized. 
Unfortunately, it is not possible, at the current state-of-the-art in man- 
machine modeling, to predict the degree of training a, pilot needs before his 
control strategy stabilizes. The given pilot data does seem to violate the 
underlying assumption in the optimal control model — that the pilot is well-train- 
ed. 


Comparing the mean sink-rate and throttle input for the pilot and model 

shows that the pilot is probably overcorrecting for C , whereas the model is 

L , §e 

undercorrecting. This results in larger pilot throttle inputs and lower sink- 
rates than the model, and a subsequent tendency for the pilot to overflare. 

(Note the rapid pilot attempt to compensate for t > 3 sec.) 

< 

Pitch angle and tail input are shown in Figure 17. The basic aircraft 
approaches the flare window in a -4° nose-down attitude. The pilot is required 
to (try to) land with a positive 2° pitch. Neither the pilot. nor the model 
is capable of achieving a nose-up attitude at touchdown. The pilot data shows 
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FIGURE 14 FLIGHT PATH DISPERSIONS, UNAUGMENTED STOL 
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FIGURE 15 SINK RATE, UNAUGMENTED STOL 
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FIGURE 16 THROTTLE/ THRUST INPUTS, UNAUGMENTED STOL 
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a much higher variability than does the model. The apparent reasons for this 
are much the same as in the autospeed case, i.e., run-to-run variability in 
pilot strategy and low attentional allocation to the pitch attitude indicator. 
The model's pitch response was found to be quite sensitive to changes in q . 

For example, q Q = 0 yields a model prediction 0 , - -2°. The mean tail input 
for both model and pilot are in close agreement with the pilot input being only 
1/4° less on the average for t > 2 sec. 

4.2.4 Sink-Rate Predictions at Touchdown 


Among the most important predictions generated by the model are the prob- 
ability density of the touchdown velocity and the mean touchdown velocity. 

These quantities are given by Eqs. (189) and (185), respectively. It should 
be re-emphasized that our predictions of touchdown statistics are based only 
on those flight paths that terminate within approximately 1 sec of the nominal 
touchdown time. The one second "cutoff" has been chosen arbitrarily; for flights 
tending to land late, model predictions are probably not highly representative 
of pilot actions. 

Figure 18 shows the model's prediction of the touchdown velocity probability 

density (conditioned on T, < 5.9 sec) for the three STOL configurations. The 

d — 

comparisons are interesting. The augmented STOL density shows a well-defined 

peak at h - -2.25 ft/sec. In the autospeed case, the density peak has shifted 

to -3.25 ft/sec indicating the higher touchdown velocities associated with 

earlier landing times (see Figure 9). In addition, the density has broadened, 

predicting more variability in touchdown velocity. The unaugmented results show 

an even greater flattening of the touchdown density, to correspond with a highly 

<A? 

variable sink-rate at touchdown. 

Because of the limited number of experimental runs that terminate prior 
to t = 5.9, it is not realistic to compute a probability density or histogram 

The sensitivity of 0 to 6t is evident by comparing the mean responses of 
Figure 17a and 17b. 

kk • • • 

The density of h increases for -h < 1.5 ft/sec since h constant - 1.35 for 
T^ > 5.4 sec. Ideally, there would be an impulse at ft = -1.35 with area = .2. 
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FIGURE 18 PROBABILITY DENSITY OF TOUCHDOWN VELOCITY, T d < 5.9 
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of h . However, it is plausible to compare the average touchdown velocity as 
computed by Eq. (183) with the average sink-rate of those data trials that land 
prior to t = 5.9. Table 1 shows this comparison. P(5.9) is the predicted 
probability of landing prior to t = 5.9. N^ is the number of data runs that 
have landed by this time. The model and data results are in reasonable agree- 
ment, despite the low N^, which gives rise to the large standard deviation in 
the data mean. 


TABLE 1 

MEAN SINK-RATE FOR TOUCHDOWN < 5.9 SEC 



MODEL 


DATA 


Condition 

P (5. 9) 

h J 
td 

N. 

l 

N. 

l 

N 

*td 

S.D. 

Augmented 

.69 

2.25 

7 

.7 

3.4 

1.7 

Autospeed 

.79 

3.10 

8 

.8 

3.3 

1.5 

Unaugmented 

.80 

3.10 

5 

.62 

5.0 

2.0 


4.3 POTENTIAL MODEL MODIFICATIONS 

When viewed in total, the results of Section 4.3 show that the model has 
the capability of predicting pilot performance in a landing task. Certain 
aspects of the model (e.g., our treatment of the terminal control aspects of 
the landing) were validated by the data. On the other hand, the model was 
found to possess certain deficiencies that indicate modifications to the model's 
a prior structure and/or parameter values. 
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4.3.1 Estimation of Ground Effect 


The data for the three conditions studied showed that as the task became 

more difficult, the pilot exhibited a greater tendency to overflare. The model, 

on the other hand, indicated a trend towards landing short. The pilot seems to 

be overcorrecting for C T while the model is undercorrecting. The under- 
age " . 

correction of the optimal control model is evident in the thrust predictions 
of Figures 7, 11 and 16. As each condition increases in difficulty, the differ- 
ences between pilot and model 5T c become more pronounced. 


The optimal control model has the ability to estimate C = x i continu- 

L* , ge J. 

ously from displayed information given an a priori estimate x^(0) = -.09. 

• • 

There is no ability to estimate x^, or to use an initial estimate of x^. 

Therefore, as the quality of the pilot’s information base decreases and as the 

• & 

task becomes more difficult, the estimate x^(t) changes less quickly. In 
terms of the model elements, the Kalman filter has larger time constants. Fig- 
ure 19 shows clearly this trend. In the augmented case, the pilot’s estimate 

of CL lags the true value with a .05 estimation error. In the unaugmented 

L.ge 

case, the estimation error has doubled; as a result, pilot compensation for 
is too small. 


From the above discussion, and the results of Section 4.2, it is apparent 

that the model should have the ability to estimate (at least) C . This is 

i»g e 

easily accomplished within the model framework by defining two states x^ = C^, 
X 2 = C L , with ... 


x^t) = x 2 (t) 

■ . I 

* 2 (t) = C L (t) = z(t) 


(194) 


The above equation (194) would then replace the equation for C in the STOL 

L 

• *• - 
model. In this case, z(t) will be a function of h(t). Since the model s 

• * , f ' 1 . r , ' 

In the limiting case of no display information x^(t), = x^(0) for all t. 
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Kalman filter generates estimates of all states x.(t), the filter will continu- 

• ^ 

ously estimate C (t). The net result would be an estimation of C that is 
t L 

closer to the actual value, but showing a tendency to overestimate since the 
Kalman fitler would always be expecting C to increase. 

J-i 

The above discussion is also appropriate to the estimation of the pitching 
moment ground effect, However, no conclusions can be drawn from the 

poor pilot pitch data, as regards the need for modifying the estimation of C^. 
Better estimation of will probably produce little effect on the model's 

pitch response due to the overwhelming effect of the STOL's poor pitch damping. 


4.3.2 Cost Functional Weightings 

The cost functional weightings for the optimal control model were chosen as 
described in Section 3.3 by considering pilot regulation about a nominal trajec- 
tory. With respect to altitude regulation, only a weighting on deviations 6h 

was included. The results of Figures 8 and 12 indicate that the pilot is probably 

• •• 

regulating about h , and possibly h , as well. Therefore, it seems 
° nom’ r J nom’ 

appropriate that any model modifications should include cost functional weight- 
• •• 

ings on Sh and (possibly) 6h. It can be shown that if these deviations were 

2 . 2 * 

weighted, it becomes unnecessary to include u = 6T in the cost functional. 

c 

The control input 6T^ is implicitly kept small by the model in order that 6h 
remain small. A reexamination of the data with emphasis on model-data matching 
or identification techniques is warranted. This is a relatively important task 
as regards further modeling efforts, since one must determine whether a pilot 
tends to control altitude deviations directly, or indirectly via control of 
sink-rate deviations. 


Changes in the cost functional weighting on 0 deviations was found to 
induce relatively large changes in the model's resulting pitch response. Since 
the model, on the average, controlled pitch somewhat better than the pilot, it 
appears that the a priori value of q Q = .25 is too large. Setting q„ = 0 

0 V 


In most past modeling efforts, it was not necessary to weight u . 
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would bring model predictions more in line with the data. However, any attempt 

to identify an accurate "value" of q Q is virtually impossible from the avail- 

o 

able pitch data. It appears that the qg associated with the pilot response is 
variable on a run- to-run basis. 

4.3.3 Human Parameter Values 


Slight changes in the human response parameters (t, p^, x^, etc.) 

generally did not have a great effect on the model predictions for the STOL 

problem studied. The only parameters for which the data-model trends indicate 

potential changes are X (or p ) and the attentional allocation parameters 

ri^ in Eq. (90). It appears that the pilot is probably devoting more attention 

to the sink-rate indicator during the flare. Our assumption that p . = l/(no. 

of instruments) is probably somewhat inaccurate. A posteriori analysis of the 

data can be used to indicate the relative values of n., since the variance 

x 

of a given display quantity will be related inversely to the associated rp . 
Methods for predicting a priori fractional allocation of attention or scanning 
behavior are beyond the scope of the immediate effort. However, this topic is 
of central importance for display design and panel layout. 

Lower values for the motor noise/signal ratio, or corresponding higher 
values for X seem to be indicated by the autospeed data. In this case, the 
model's control variability is greater than the pilot's. It should be empha- 
sized that there is relatively little data available in the manual control field 
pertinent to multiple control input tasks. We have assumed values for X - .1 
on the basis of past single-axis task studies. In multiple control situations, 
it may be reasonable to expect increased values of x depending, for example, 
on the relative use the pilot makes of each control.* 


The concept of pilot control allocation is interesting, and may have a counter- 
part to the display allocation problem. 
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4.3.4 System Modeling Errors 


In developing the human operator model, it was necessary to derive a 
linearized state-space model for the vehicle dynamics. These dynamics are 
highly nonlinear with respect to angle-of-attack variations, thrust variations, 
etc. Thus, slight differences between the model predictions and pilot data 
might be attributed to imprecise linearization. For example, changes in the 
pitch damping terms in the linearized model have a pronounced effect on the 
short period mode; changes in the thrust characteristics will change the manner 
in which in which the pilot applies control, etc. In addition to the above, 
there are modeling errors introduced in our simple modeling of the autospeed 
and pitch command-and-hold systems. When the autospeed system is engaged, we 
have set u' = 0; the pitch hold system implied 66 = 6q = 0 in our model. 

In actuality, these augmentation systems are not perfect and introduce changes 
in the model states. For example, the pitch hold system in the fully augmented 
case was found capable of controlling 8 to within only .5° during flare. A 
more extensive model analysis could include the dynamics of the augmentation 
systems, but the price to pay is an increased model complexity. 



V. CONCLUSIONS 


Optimal control and estimation techniques have been combined with human 
response theory to develop a model for the pilot in a STOL landing task. This 
model extends previous work in man^machine systems analysis where only the 
approach portion of the landing task was considered. The ability to carry model 
predictions through flare, to touchdown, now exists. 

i; - . 'I’ 

The model development began with the nonlinear, longitudinal equations of 
motion for the aircraft. These equations w;ere linearized about the equilibrium 
flight conditions, to develop a linear model for the STOL dynamics suitable for 
analytic investigation. The linear model coefficients were obtained from NASA 
wind tunnel measurements. The ground effects, which are an important aspect of 
the STOL landing problem, were included in the model equations. 

The optimal control approach to human operator modeling provided the basis 
for developing the pilot model. The underlying assumption in this approach is 

‘ ' V •* * 

that the well-trained pilot will behave in an optimal manner subject to his 
inherent limitations and the task requirement. Thus, the model that has been 
developed includes such pilot limitations as time-delay, inherent randomness, 
visual threshold effects, display attentional allocation, and "neuro-motor" 
dynamics. The terminal control aspects of the landing task were modeled by 
assuming that the pilot regulates the aircraft motion about a nominal flight 
path that meets the desired touchdown conditions. The nominal path has the 
property of being generated with minimal control effort. 

The pilot model is time-varying, with the time-variation arising from the 
terminal control aspects of the task and the pilot's estimation of, and adapta- 
tion to, the ground effects. The model's output predictions are of a covariance 
propagation type, i.e., for any system variable, we obtain the mean and standard 
deviation as a function of time from flare initiation. The mean and standard 
deviation are sufficient statistics for predicting the joint probability density 
of the system state at any given time. These probability densities were trans- 
lated into probability densities at touchdown to give model predictions of 
touchdown time and touchdown velocity. 
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The model was applies to predict flare performance for three different 
levels of STOL augmentation: 

1. autospeed plus pitch' command and hold, 

2. autospeed system only, and 

-?■ 3. no augmentation. 

Parameter values were chosen a priori on the basis of the task requirements and 
existing results in human response theory. The same set of model parameter 
values were used across the three STOL cases studied. The model predictions 
were compared with pilot data obtained on the NASA-LRC Real-Time Dynamic Simu- 
lator. The comparisons were generally in agreement, with the most accurate 
predictions being obtained for the fully augmented aircraft. The data indicated 
that certain aspects of the model were valid, e.gl, our method for treating the 
terminal control problem, but that other facets require slight modification, 
e.g. , the pilot's estimation of the ground effect. Some suggested modifications 
to the model, that will enhance its predictive capability, were discussed. 

The pilot landing model can serve as a very usdful design tool in analytic 
investigation of landing performance. It is relatively easy to study the effects 
of changing the display information, display gains, etc., thus obtaining prob- 
abilistic predictions associated with different display formats. Adding, or 
deleting, vehicle stability augmentation systems can be studied merely by chang- 
ing the vehicle parameters; the parameters associated with the pilot (.e.g., 
time-delay, cost functional weightings, etc.) remain the same. Thus, it is 
possible to exercise the computer model in a straightforward, repetitive manner 
to study' the effects on performance of vehicle and/or display modifications. 
Extensive simulation programs can be eliminated, since the model can serve to 
glean the most promising configurations for subsequent pilot simulation tests. 

Further efforts in extending the modeling techniques that have been devel- 
oped include the following: 
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1. Generation of Computer Simulation Flight Trajectories 

The predictions generated by the optimal control model are of a prob- 
abilistic nature. Thus, with a single computer - run, we generate the 
statistics associated with an entire ensemble of trajectories. However, 
in certain applications it may be desirable to obtain representative 
sample flight paths in order to investigate potential oscillations, 
overshoots, etc., that might not be readily apparent from the averaged 
statistics. The ability to use the model in a "Monte Carlo" manner 
would certainly enhance its flexibility, and is a relatively straight- 
forward extension to the present effort. 

2. Display Attentional Allocation 

As discussed earlier, it is important to include within the pilot model 

the capability to predict a pilot's fractional attention allocation 

among a set of display indicators. This is a necessary step towards 

the eventual use of the pilot model for display design. One method 

of approach is to assume that pilot scanning (or attention) strategy 

is chosen to minimize the uncertainty in estimating those state 

T31 

variables most important for feedback control. Accurate pilot 

scanning data will be necessary for validating any such scanning sub- 
model that is eventually developed. 

3. Pilot Control in Terminal Areas 


The model that has been developed is appropriate for the specific 
pilot task of landing approach and flare. However, the methodology 
we have used can be extended to model a pilot's ability in following 
curved approach paths (or multi-segment paths) appropriate in a 
terminal configured vehicle context. Thus, one would have the ability 
to determine the degree of pilot error in executing given maneuvers, 
including compensation for changes in relative wind, timed arrival 
at given waypoints, etc. 
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4. Optimization of Display Information 


One of the ultimate goals of the pilot model is its application to 
optimizing the information content of a display panel. To achieve 
this goal, it will he necessary to develop (either via experiment or 
analytically) more precise relationships between a display symbol and 
related model parameters such as observation noise, perceptual thresh- 
olds, etc. A' model can then be used to rank-order proposed display 
panel formats, with model predicted performance serving as a figure 
of merit. 

5. Display Monitoring 

The. model's, information base consists presently of those instruments 
used by, the pilot in generating a given control input. There is no 
explicit- treatment of pilot monitoring displays that are not directly 
appropriate for control, but may relate to vehicle "status". It will 
prove important to model these peripheral sources of information when 
considering pilot identification of • system failures, anc the broader 
aspect of total flight management. The model has the capability to 
treat such "status" information as discrete information, obtained 
periodically, as opposed to "control" information that is processed 
continuously. Modifying the model to include continuous and discrete- 
time information is an area for further effort, especially as regards 
to pilot monitoring in automated and/or semi-automated systems. 
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APPENDIX A 


SOLUTION OF THE LINEAR OPTIMAL TERMINAL CONTROL PROBLEM 


This appendix derives the solution to the linear regulator problem where 
it is required to meet prescribed set of terminal conditions. 

Problem Formulation 


The system dynamics are described by . the linear time -invariant equations 

x(t) = Ax(t) + Bu(t) (A.l) 

x(t ) = x = given (A. 2) 

o o 

where x = col [x, ,^ 2 * • • • » x ]• is the system state and u = col t u i> u 2» • * ,u r ] 
is the control. It is desired to drive the sytem from the initial state x 

o 

to a terminal state x(t^) that satisfies the p boundary conditions 

Hx(tj) + c = 0 (A. 3) 

where H is a given p x n matrix of full rank p £ n and c is a p-vector. 
In addition, the control that accomplishes the transfer must minimize the cost 
functional 


J(u,x 

o 



(t^)Fx(t 



[x' (t)Qx(t) +u' (t)Ru(t) ]dt 


t 

o 


(A. 4) 


where F and Q are positive semi-definite matrices and R is positive 
definite. 


It is assumed that there exists at least one control that can drive the 
system from x q to an x(t^) satisfying (A. 3) so that the optimization problem 
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posed is non-vacuous. If the system (A. 1) is completely controllable, this 
assumption will be satisfied. 

Problem Solution t - • > i 

i \ 

In order to solve the optimization problem, the constraint equation (A. 3) 

is appended to the cost function J(u,x ,t ) by introducing a Lagrange multiplier 

o o 

V. Thus, we define a new cost functional, 

J(u,x Q ,t o ) = (Hx(t^) + c)'v + J(u,x o ,t Q ) (A.5) 

‘ , 4 V* ' 

where V will be chosen such that the terminal condition Hx(t^) + c will be 
small — hopefully zero. 

The minimization of J subject to Eq. (A.l) is easily accomplished using 
Hamilton Jacobi theory. ^ Define the Hamiltonian 

H^x, u) = x’Qx + u’Ru + -|^ (Ax + Bu) (A. 6) 


where the function J(x,t) of x and t is the optimal cost 

J(x,t) = min J(u,x,t) (A. 7) 

u 

for the system (A.l) with initial state x at time t. Clearly, at t = t , 

J(x,t^) = x'Fx + (Hx + c)'v (A. 8) 

Jc 

The optimal control, u , must minimize the Hamiltonian. Thus, 

|S . o (A. 9) 
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which yields 


The Hamilton- Jacobi theory then states that the function J(x,t) satisfies 
partial differential equation 

M + -(«. I- •*) - ° 

with the boundary condition (A. 8). Substituting Eq. (A. 10) in (A. 11) gives 


3J 3J' 
3t 3x 


The solution J(x,t) 


Ax + x'Qx - j |^' BR -1 B' = 0 , 

can be shown to be, by direct substitution. 


J(x,t) = x'K(t)x + x'h(t) + r (t) 
where the matrix K(t), t £ t^ satisfies the Riccati equation 
K(t) = - K(t)A + A'K(t) + K(t)BR -1 B' K(t) - Q 


with boundary condition 
K(t x ) = F 

Note that if F = Q = 0, then K(t) = 0. 

The function h(t), t < t^ is generated by 

h(t) = - (A - BR _1 B' K(t) ) * h(t) = - A’ (t) h(t) 


Mtp = H'V 


(A. 10) 
the 

(A. 11) 

(A. 12) 

(A. 13) 
(A. 14) 
(A. 15) 

(A. 16) 
(A. 17) 
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Finally, r(t) = v’c. Hence, the optimal control is 


u(x,t) » - R~ 1 B , K(t) x (t) - R~ 1 B , h(t) 


It remains to determine v such that the boundary condition (A. 3) is 
fied. We note that the solution for h(t) is given by 

h(t) = (t 1 ,t)H'v 

where $(t,r) is the state transition matrix corresponding to A(t), i.e. 


— <Kt, t) = A (t) $(t, t), $(t,T) = I + 


The state trajectory of Eq. (A.l) is then generated according to 


x(t) = (A - BR 1 B'K(t)) x(t) - BR VhCt) 


- A(t) x(t) - BR _1 B' $'(t 1 ,t)H , V 


x(t^) is then given by 


x(t n ) = $(t ,t )x - W(t,,t )H'v 
i loo 1 o 


where 


W(t ,t ) = 
l o 


r 


$(t lf T) 


BR 1 B' 


(t lf T) dx 


+ If K(t) = 0, $(t,x) = e A(t " T) . 


(A. 18) 
satis- 


(A. 19) 


(A. 20) 


(A. 21) 


(A. 22) 


(A. 23) 
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The matrix W(t- ,t ) 
lable.' 7 ’ 


is invertible if the system (A.l) is completely control- 


Premultiplying Eq. (A. 22) by H it can be seen that for Hx(t^) = -c, 
the vector V should be chosen as 


V = [H W(t 1 ,t o )H']' 1 * [H $(t 1 ,t o )x o + c] 


(A. 24) 


The required inverse exists if the system is completely controllable, and since 

£ 

H was assumed to have full rank. This completes the solution to the originally 
posed problem. 

Extension to Control Rate Weighting 

A 

In the optimal control of human response, the cost functional J(u,X Q ,t o ), 
Eq. (A. 4), is of the (general) form 


J(u,x o ,t Q ) = x' (t 1 ) Fx(t^) 


+ 



(x'Qx + u’Ru + u'Gu)dt 


(A. 25) 


The system dynamics and boundary conditions (A.l) - (A. 3) remain unchanged. In 
order to apply the above results to solve this problem, we define an augmented 
n + r state vector y = col [x,u] and define y(t) = u(t). Thus, the equation 
(A.l) becomes 

y(t) = A Q y(t) + B Q y(t) (A. 26) 


* ( V ' *o 


(A. 27) 


The assumption of full rank is tantamount to having non-conflicting boundary 
conditions. 
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where 


'A { B~ 


' 0 ■ 


r x -) 

i 

1 

> B o = 



0 


’ X o = 


i — 

o 

o 


I 

-J 


0 


(A. 28) 


The cost functional Eq. (A. 25) and terminal condition Eq. (A. 3) may then be 
written as 

rh 

J(y,X 0 .to> = x(t x ) + / (x'Q 0 + y'Gy)dt (A.29) 

t 

o 


HqXC^) + c = 0 


where 


(A. 30) 


F q = diag (F,0) , Q q = diag (Q,R), H q = [H,0], (A.31) 

The derived results of Eqs. (A. 14) - (A. 24) are now applicable to the prob- 
lem of minimizing control rate, with the replacements A ■+ A q , B B q , Q •+■ Q q , 

R -*• G, H-*H , x X> u' y y = u, etc. 
o 


no 



APPENDIX B 


NUMERICAL SPECIFICATION : OF THE STOL MODELS 


The elements of the state space STOL models are functions of velocity, 
physical characteristics of the aircraft, and the aerodynamic coefficients. These 
elements are analytically specified by Eqs. (57) - (61) in the text. In this 
appendix, the numerical values for the A and _B matrices are presented for the 
basic aircraft, the aircraft with autospeed, and the aircraft with pitch command 
and autospeed. 

Basic Aircraft 


The basic aircraft approaches at an initial airspeed of 75 knots on a 6° 
glideslope. The centerline of the aircraft is initially pitched down -4.1°. The 
A and B matrices are given by: 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-.0373 

.0855 

0.0 

-.255 

0.0 

. 0.0 

-.063 

0.0 

-.522 

-.566 

1.0 

.0268 

0.0 

-.134 

0.0 

.63 

.2087 

-1.2 

-1.47 

-.0107 

0.0 

-.0212 

0.0 

0.0 

0.0 

0.0 

1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-.105 

-1.0 

0.0 

1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-2.0 
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0.0 

0.0 

0.0 

0.0 

0.0 

-.0043 

0.0 

-.096 

0.0 

-3.5 

0.0 

0.0 

0.0 

0.0 

2.0 

0.0 


Aircraft With Autospeed 

The aircraft with autospeed approaches at an airspeed of 75 knots on a 6° 
glideslope. The airspeed is held at 75 knots and the centerline of the aircraft 
is initially pitched up 2°. The A and B matrices are given as 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-.063 

0.0 

-.556 

1.0 

.0268 

0.0 

-.134 

0.0 

0.0 

-1.2 

-1.47 

-.0107 

0.0 

-.0212 

0.0 

.63 

0.0 

1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-1.0 

0.0 

1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-2.0 
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i 

The fully augmented aircraft approaches on the 6° glidelsope. The airspeed 
is held at 75 knots and the pitch is held at 2°. The A and 15 matrices are 
given as 



In all cases, the aircraft mass is 1700 slugs and the moment of inertia 

5 2 

about the y-axis is 2.46 x 10 slug-ft . 
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